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In all cases, agreement between theory and experiment is encouraging. 
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SUMMARY 


A frequency-domain, noncompact-source theory for both the steady-loading and volume- 
displacement (thickness) noise of high speed propellers has been developed and programmed. 
Both near-field and far-field effects have been considered. This code utilizes blade surface 
pressure distributions obtained from a 3D (three-dimensional) nonlinear aerodynamic flow 
field analysis program as input for evaluating the steady-loading noise. 

In order to investigate the effects of the aircraft installation on the sound field generated 
by the propeller, simp lified mathematical models of the velocity fields induced at the propeller 
disk by nearby wing and fuselage surfaces and by angle-of-attack operation have been 
developed to provide estimates of the unsteady loading imposed on the propeller by these 
potential field type interactions. These unsteady blade loadings have been coupled to a 
chordwise compact propeller unsteady-loading noise model to provide predictions of 
unsteady-loading noise caused by these installation effects. 

Both the steady loading and thickness model and the unsteady loading model described 
above have been checked ag ains t NASA-supplied data for the SR-2 and SR-3 model propfans. 
In addition, the steady loading and thickness model has been compared with data from the 
SR-6 model propfan. These theoretical models have been employed in the evaluation of the 
SR-7-powered Gulfstream aircraft in terms of noise characteristics at representative takeoff, 
cruise, and approach operating conditions. 

Finally, an analysis to estimate the corrections to be applied to the ffee-field noise predic- 
tions in order to arrive at the measurable fuselage sound pressure levels has been formulated 
and progr amm ed. This analysis considers the effects of fuselage surface reflection and diffrac- 
tion, together with surface boundary layer refraction. 

In all cases, agreement between theory and experiment is encouraging. 
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NOMENCLATURE (Section 2.1) 


ao Speed of sound 

B Number of blades 

Bd Chord: diameter ratio 

Cl Blade section lift coefficient 

Cp Power coefficient 

Cu Absolute tangential velocity (swirl) 

c Blade chord 

Co Speed of sound 

D, Dt Diameter 

FA Face alignment (defined in Figure 1) 
f Force 

ft Normalized chordwise blade loading function (defined in Equation 2. 1.4) 

H Normalized chordwise blade thickness distribution (defined in Equation 2. 1.4) 

i V3 

J Advance ratio of propeller 

J( ) Bessel function 

K Arbitrary constant 

kx/ky Dimensionless chordwise wave numbers (defined in Equation 2.1.3) 

L' Blade section lift 

Mht Helical tip Mach number 

Mr Relative Mach number 

Mt Rotational tip Mach number 

Mo, M x Axial Mach number 

MCA Midchord alignment (defined in Figure 1) 

m Blade passing frequency harmonic number 

n Rotational speed in rps (radians per second) 

P Power 

P LM Steady loading (lift force) source 
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Pv M Volume displacement (thickness) source 

p Pressure 

R/r Blade radius or observer distance (defined in text) 

r 0 Observer distance 

T' Torque per unit span 

tb Maximum thickness: chord ratio 

u/v/w Velocities 

W Power 

X Nondimensional chordwise coordinate, equals x/c 

x Chordwise or axial coordinate 

y Sideline distance 

z Radius ratio (r/rtip) 

3 P . Blade pitch angle (defined in Figure 1) 

7 Ratio of specific heats 

A Differential 

e Emission angle 

Ol Defined in Sketch 2 . 1 

d Swirl coefficient (Cu/rn) 

po Density 

4> Phase shift (defined in Equation 2. 1 . 5 ) 

4* Blade thickness and loading distributions in the frequency domain (defined in 

Equation 2 . 1 . 4 ) 

n Rotational speed in radians per second 
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NOMENCLATURE (Section 2.2) 


B Number of propeller blades 

c Ambient speed of sound 

fzo Steady (thrust) force per unit spanwise length exerted by each blade in the z direc- 

tion (Figure 37); fzo is a function of radius “r” 

f(j)o Steady (tangential) force per unit spanwise length exerted by each blade in the “<j>” 
direction (Figure 33); note that f^ois a function of radius and also that, in view of 
Figure 33, if the sense of rotor rotation is clockwise (forward looking aft), f<j,o will 
be negative 

h Representative airfoil thickness as a function of radius 

j VT 

J Bessel function of the first kind 

1 Airfoil chord as a function of radius 

Mf Flight Mach number 

n Harmonic of blade passing frequency of interest (n s 1) 

p 0 Far field acoustic pressure due to steady loading and thickness 

r Radius of interest 

R Distance of observer from center of retarded propeller disk location (reference is 

made to Figure 38) 

Rc Distance of observer from center of current propeller disk location (Figure 38) 

t Time 

$ Azimuthal angle as defined in Figure 33 (note that 4> = 0, 180 represent the 

horizontal direction) 

k Wave number, nB 0/c 

po Ambient air density 

e “Emission” angle for propeller noise 

0c “Observer” angle for propeller noise 

Cl Angular velocity of propeller in radians per second 
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1.0 INTRODUCTION 


1.1 Background 

Recent studies have shown that the high speed turboprop offers an attractive alternative 
to turbofan propulsion for certain future commercial aircraft applications because of its sub- 
stantially higher propulsive efficiency potential. However, one of the major obstacles to the 
use of high speed turboprops is the high near-field noise level generated by the propeller, 
which may cause a potentially serious cabin noise problem. In response to this potential 
problem, there has been substantial activity in both theoretical and experimental investiga- 
tions of the noise characteristics of high speed propellers by NASA and private industry. 

The currently envisioned high speed turboprop systems call for cruise flight speed Mach 
numbers of 0.7 to 0.8 and propeller tip speed Mach numbers of 0.7 to 0.9, which results in 
propeller blade operation in the transonic and supersonic regimes. Recent theoretical models 
of propeller noise have been produced to account for forward flight effects, source non- 
compactness, and supersonic motion of the blade-fixed sources relative to the observer 
position. These effects are all present to some degree for propellers operating at transonic and 
supersonic relative Mach numbers. The concept of the propfan, a low aspect ratio, high blade 
number propeller with swept blades, has features which promise significantly lower noise than 
conventional propellers. These features include very thin blade sections with large sweep 
angles near the tip to produce subsonic relative Mach number components normal to the lead- 
ing edge of the blade, thus reducing losses and flow field discontinuities associated with the 
formation of shock waves. 

Experimental measurements of several propfan configurations under NASA (and NASA- 
sponsored) programs have shown that propellers operating at transonic tip relative Mach num- 
bers in actual (and simulated) forward flight exhibit noise characteristics considerably different 
from those predicted by earlier propeller noise theories. These differences include the trends 
in level with increasing tip relative or helical Mach, number, as well as the harmonic spectrum 
shape. Since the newer theories predict the measured trends quite well in the subsonic regime, 
there is reason to suspect the failure in the theories to be related to the changes in character 
of the flow in the transonic regime being not adequately simulated by the theories as they now 
stand. 

The substantial improvements in aircraft fuel efficiency projected for turboprop systems 
over turbofan systems gives strong incentive to resolve the discrepancy between predicted and 
measured noise trends in the transonic/supersonic operating range. If an accurate prediction 
model which contains the essential physics of the noise production processes involved can be 
established with confidence, then the procedure can be exploited to understand the noise 
production mechanisms and their controlling parameters. From this understanding, a design 
optimization procedure can be developed which yields the best combination of high efficiency 
and low noise. An accurate prediction of the noise characteristics of a propeller, especially its 
near-field spectral and directivity characteristics, will also permit a better chance to design the 
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most effective cabin structural and airborne noise attenuation features and, hence, m i n i m ize 
cabin interior discomfort with the smallest possible impact on the weight of the aircraft. 

Although the issue has not been resolved, several reasons and/or explanations for the 
apparent failure of current theories of propeller noise in the transonic/supersonic regime have 
been put forth. One possible explanation (which we believe to be the correct one) is that the 
aerodynamic input to the acoustic theory has thus far been inadequate. Current practice has 
been to input blade surface pressure distributions using a polynomial distribution matched to 
an overall lift force at each spanwise location along the blade, or to utilize linearized two- 
dimensional strip theory analysis to calculate the input pressure distributions. Current noise 
theories utilize this input pressure distribution to compute the contribution to the noise field 
due to blade loading. The blade geometry, blade cross section shape, or thickness distribution, 
is used to compute the blade volume-displacement or thickness contribution to the noise field. 
Both loading and thickness noise produce tones at multiples of blade-passing frequency. The 
thickness noise component is a function of the spatial/time variation of the component of blade 
surface velocity normal to the surface itself. This is a function only of blade shape and blade 
motion relative to the observer, and so is independent of the aerodynamic flow field around 
the blades. Hence, if the thickness noise component of the theory is in error, it is due to an 
inadequacy of the acoustic theory itself, and not due to an inadequate aerodynamic input. 

The loading noise component of the theory, however, could possibly be in error as a result 
of shortcomings in the aerodynamic input The inadequacy can be associated with at least two 
major effects: 

1. Three-dimensional nature of the flow in the tip region of the blade, characterized 
by spanwise flow migration and blade tip “relief* effects 

2. Nonlinearity of the flow in the transonic regime, including the effects of shock 
formations and supersonic “patches” on the blades. 

The first of these effects is present even in subsonic flow, where the current acoustic 
theories seem to work well enough. The second effect is unique to the flow regime of interest, 
and it is strongly suggested as a necessary effect which must be properly modeled in order to 
predict the noise characteristics of propellers in the transonic regime properly. Although the 
three-dimensional effects may alter the nonlinear and shock-related features of the blade flow 
field, it is strongly felt that it is the nonlinear effect which may resolve the data/theory 
discrepancies in the transonic regime. Although three-dimensionality is certainly present, its 
effect by itself is probably not important since the two-dimensional aerodynamic input 
provides reasonable results in the subsonic regime. 

It can be argued that future turboprop propellers will be of the propfan type, with swept 
blades of thin cross section, and hence will be “designed” for minimal transonic nonlinearity 
effects and shock-free flow. However, even if this situation could be achieved at “design” 
conditions, the aircraft propulsion system is required to operate at off-design conditions during 
takeoff, climb, descent and landing approach modes, and thin blade sections operating even 
slightly off-design will exhibit nonlinear flow distributions with shock formations. These 
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off-design conditions may well constitute a major fraction of the flight duration for the typical 
500-mile missions (or less) envisioned for high speed turboprop application. 

There are other possible causes of the disparity between current theories of propeller noise 
and recent measurements. These include: 

• Nonlinear propagation effects due to the impulsive nature of the blade-generated 
wave forms in the transonic regime 

• Other sources such as nonuniform loading noise due to installation-imposed distor- 
tions on the propeller 

• Fuselage boundary layer refraction of the incident sound 

• Shielding of the incident sound by nearby wing surfaces. 

Some of these may not apply to all the recent experimental observations acquired by NASA 
but, nonetheless, are effects which need to be considered in evaluating a given turboprop 
aircraft installation. 

1.2 Objectives 

The primary objective of this program has been to provide an improved, unified theoretical 
model for predicting the noise of high speed propellers which can satisfactorily explain 
observed experimental trends in the transonic helical Mach number regime (as discussed in 
Section 1.1). A second objective of this program was to provide preliminary theoretical models 
for predicting the installation effects on propeller noise. These effects include additional 
unsteady loading noise due to angle-of-attack, wing and pylon/nacelle interference effects; and 
aircraft surface reflection, refraction, and shielding. The primary objective focused on an 
improved understanding of the basic propeller noise generation due to steady loading and 
volume-displacement (thickness), while the secondary objective focused on additional, 
installation-induced noise sources and modification of the basic noise source characteristics 
by the installation (aircraft) environment. 



2.0 ANALYSIS AND RESULTS 


2.1 Steady Loading Theoretical Model 

2.1.1 Introduction 

2.1 .1.1 Acoustic 

The starting point of the acoustic analysis is described by the so-called “Ffowcs Willi ams- 
Hawkings Equation” (Reference 2.1.1): 


4np (x,t) * 3 \ j r|l-M r l dS ‘ 9^ f\_ r|l-M r l] 


dS 


( 2 . 1 . 1 ) 


where the square brackets imply evaluation of their contents at the “source” or “retarded” 
time t, given by 



when t = “observer” time. It should be noted that Equation 2.1.1 is quoted omitting the 
aerodynamic stress or quadrupoie terms. 

It can be seen that if Equation 2.1.1 is evaluated as is, namely in the time domain; singu- 
larities will occur when M r , the Mach number of the source in the direction of the observer, 
becomes unity. For this reason, many workers (for example. References 2.1.2 through 2.1.5) 
have taken the step of Fourier Transforming Equation 2.1.1 into the frequency domain in order 
to make use of this formulation in the prediction of the sound generated by an open rotor. 

Hawkings and Lowson (References 2.12 and 2.1.5) employed Equation 2.1.1 in the calcu- 
lation of sound from open rotors operating at supersonic tip speeds. They assumed that the 
blades were thin and that the surface integrals could be replaced by integrating over the mean 
planform area. The thin blade assumption also led to the dropping of the quadrupoie term 
from the original equation. 

Hawkings and Lowson compared their predictions with data measured by Hubbard and 
Lassiter (Reference 2.1.6) and Kurbjun (Reference 2.1.7) and concluded that discrepancies 
in the comparisons were the result of nonlinear propagation of the acoustic signal. They 
employed the weak shock theory of Whitham (Reference 2. 1.8) to account for this nonlinearity 
and were, by this means, able to improve the agreement between the reported data and the 
predicted results by a considerable degree. It is, perhaps, worth noting here that the data were 
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taken in the plane of rotation of a two- and a three-bladed propeller at distance-to-diameter 
ratios of 7.66 and 10.00, respectively. In addition, the information available concerning blade 
section geometry and loading was sparse. Several points made by Hawkings and Lowson do, 
however, bear repetition here. They observed that at supersonic source speeds it is possible 
for the source to emit signals from numerous locations which may arrive at the observer simul- 
taneously. They also noted that changes in the lift or thickness distributions across the blade 
chord have direct and predictable effects upon the acoustic harmonics and that this suggests 
the possibility of altering the blade profile and aerodynamics to produce an acceptable acoustic 
spectrum. 

Hawkings and Lowson used a cylindrical polar coordinate system in their derivation. They 
also neglected the effects of forward flight - a reasonable simplification since both sets of data 
used for their comparisons came from propellers on test stands. In their laterpaper (Reference 
2.1.5) they do perform the extension to forward hub motion, observing that, acoustically, the 
effect of the hub velocity is equivalent to increasing the rotational speed by the Doppler factor 
(1 - Mo cos 0)* 1 while simultaneously reducing the observation distance to r (1 - M 0 cos e). 

The above-discussed approach of Hawkings and Lowson can be traced back to the Lighthill 
theory of aerodynamic sound (Reference 2.1.9) by means of the work of Ffowcs Williams and 
Hawkings (Reference 2.1.1). 

Taking the same starting point, but following Goldstein (Reference 2.1.10) in the use of 
generalized Green’s Function solutions, Hanson (Reference 2.1.11) defined a space-fixed, 
locally orthogonal coordinate system, tied to the helicoidal surface swept by the blade pitch 
change axis as the blade rotates and advances. The advantage of this coordinate system lies in 
the ease with which conventional airfoil coordinates can be input into the computational 
scheme. 

At the start of this work it was decided to develop and program a linearized, distributed- 
source, far field, frequency-domain model for steady loading and thickness noise only (fol- 
lowing the arguments of Hawkings and Lowson), in which the helicoidal surface coordinate 
system of Hanson would be employed. It was decided that, following comparisons with data 
taken in the acoustic near field, extensions to the theory would be formulated to account for 
near field effects. 

2.1. 1.2 Aerodynamic 

Use of a distributed-source prediction procedure requires knowledge of both spanwise 
and chordwise distributions of blade loading and thickness. The propfan designs under consid- 
eration for application on high speed commercial transport aircraft consist of rotors containing 
many highly-swept blades, in stark contrast to the two- and four-bladed propellers that can be 
seen on today’s general aviation aircraft. In addition, the aspect ratio (span/chord) of the prop- 
fan blade is small by comparison with a conventional propeller. These features combine to 
produce high values of solidity (chord/spacing) for the blading, especially in the lower portions 
of the span. In these regions, blade-to-blade interference effects cannot be neglected, and in 
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effect, the flow field in this region bears closer resemblance to that in a turbomachine than to 
that experienced by a conventional propeller. 

Advances in computational fluid dynamics and in high speed computers in general have 
led to the development of sophisticated computer codes for solving the nonviscous Euler 
equations of fluid flow in three dimensions. The first of these to gain general acceptance in 
the turbomachinery world was written by Denton (Reference 2.1.12). Since the appearance 
of the original Denton code, many other solvers have appeared (for example. References 
2.1.13 and 2.1.14), each offering subtle differences in the assumptions made, the grid 
employed, and the computational scheme involved. Since viscosity is the only flow 
phenomenon neglected in these codes, nonlinearities in the flow (such as shock waves) are 
present in the solution. It was felt, at the start of the work described here, that the combina- 
tion of a linear acoustic theory with blade loads obtained from a nonlinear aerodynamic solu- 
tion would provide an excellent basis for prediction of the sound generated by a high speed 
propfan of the type under consideration for high speed transport applications. 

In addition, it was considered that results from a 3D (three-dimensional) flow solution 
could be used to obtain the fluctuating pressure on a cylindrical surface surrounding the rotor. 
This pressure could then be Fourier Transformed to provide near field noise levels. The main 
disadvantage to using this technique was believed to be the amount of detail required to define 
the flow field adequately. Fine detail requires a large number of grid points, with a 
corresponding increase in computation time. 

2.1.2 Model Development 
2.1 .2.1 Far Reid Acoustic Model 

Following the arguments outlined above, it was decided that the helicoidal coordinate 
system developed by Hanson (Reference 2.1.11) would allow the easier input of blade 
geometry in standard coordinates. Consequently, the following equation was taken as the 
starting point for the program development (Reference 2.1.11, Equation 36). 



6 



In this equation: 

Pvm = volume displacement or thickness source, 

PLm = steady loading (lift force) source, 
e * emission angle 

Hd = fl/(l - M x cos e) where ft = rotational frequency in radians per second 
r = observation distance 

kx and ky are dimensionless chordwise wave numbers given by: 


2mB B D Hj. 
” M (1 - cos 0) 


and 


k 


y 


2mfi Bjj 

z M 

r 


( M r 2 cos0 - M \ 
1 - M x cos0 J 


(2.1.3) 


Bd = chord/diameter ratio. 

The parameters <l»v (k x ) and 'J'l (k x ) define the chordwise thickness and loading distribu- 
tions in the frequency domain. They are given by: 


and 


<V * 


J H (X) e x dX 


r 


/ ik X 

f L (X) e x dX 


( 2 . 1 . 4 ) 
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and are Fourier Transforms of the chordwise normalized thickness and loading distributions, 
respectively. 

The phase shifts, <{* and <jx>, resulting from sweep and offset respectively, are given by: 


^ B MCA 

= M (1 - M cos0) D 



MCA and FA or midchord alignment and face alignment, are defined in Figure 1. They 
describe the deviation of the blade section midchord from the helicoid swept out by the pitch 
change axis. 

It may be worth observing here that Hanson’s observation; namely, that the Hawkings and 
Lowson theory can be recovered exactly from Equation 2.1.2 by setting the flight Mach 
number, Mx, to zero; should be treated with caution. The factor k y above is used to split the 
lift on the blade into the thrust and drag components employed by Hawkings and Lowson and, 
essentially, uses the helicoidal vector diagram (Figure 1) to define the sine and cosine of the 
angle pL. This relation should be preserved, and thus, the decomposition of the lift force into 
its components should be carried out before Mx is set to zero. 

An alternative method of establishing the relationship between the Hawkings and Lowson 
theory and Equation 2.1.2 is to follow the approach suggested by Hawkings and Lowson in a 
later paper (Reference 2.1.5). They observed that the effect of forward flight on the acoustics 
of an open rotor is equivalent to deriving the equation for a stationary source, then increasing 
the rotational speed by the factor (1 - Mxcos e)* 1 while simultaneously reducing the observation 
distance (r) to r (1 - Mx cos e). It should be noted that the reduction in observation distance 
is applied to the amplitude terms only. 

In the computer code delivered under this Contract, an alternative form for the chordwise 
wave number k y has been included as an option. In this formulation: 

2mBB D M T rt (1 - M x cos0) 

= , — ■— cos p„ cos0 - — 

y (1 - M x cos0) V p zMt 



where pp is the blade pitch angle and is defined in Figure 1. 
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FA * Face Alignment (Offset or Dihedral) 

MCA * Midchord Alignment (Sweep) 

Suction (Camber Side) Surface and Pressure 
(Face Side) Surface Coordinates (y'/c)ss and 
(y'/c)ps Specified as a Function of X * x'/c 
from X ■ -0.5 to X * +0.5 


Figure 1. Definition of MCA and FA. 
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This expression uses the blade angle, 3P, rather than the helicoidal angle, 0L, to resolve 
the lift on the blade section into its thrust and drag components and, generally, has been found 
to give better agreement between prediction and data than is shown by the use of Equation 
2.13. 


2. 1.2.2 Aerodynamic Model 

In order to use Equation 2.13 to predict the tone noise generated by an open rotor, it is 
necessary to define both the thickness distribution of the blading and also the lift on the blades 
in both the chordwise and spanwise directions. 

At the outset of this work, it was decided that, in view of the fact that high speed propfan 
blading will experience relative Mach numbers in excess of unity leading to noniinearities such 
as shock waves in the flow, loading distributions generated by a 3D Euler equation flow solver 
should be used. 

The Euler equation (the Navier-Stokes equations, minus the viscous terms) can be written 
as: 


P 



3u . du . 

p ir*' >u k3i; 


| f. 

dx. x 


( 2 . 1 . 6 ) 


The general approach for solving these equations in turbomachinery applications has been 
to follow the time-marching technique of Denton. Since Denton's original work, many 
alternative computational schemes have been devised, but they are all designed for the basic 
problem of axial turbomachinery; the flow is confined within an annular duct. In order to 
extend these schemes to the case of highly swept propfan blading, two options are available. 
First, an existing code may be adapted by removal of the casing boundary to a radius much 
larger than that of the blade tip, to simulate an open rotor. Alternatively, the propfan analysis 
may be considered as a problem in its own right, and the calculation procedure and computa- 
tional grid developed to deal exclusively with open rotor geometries. 

For the purposes of this study, one example of each type of code was available. The first, 
developed at GE-CRD (GE Corporate Research and Development) is based on a code used 
extensively within GE for turbomachinery applications. It employs a Jameson computational 
algorithm (Reference 2.1.14) and is described in Reference 2.1.15. The second code, 
designated NASPROP-E, was developed under contract from NASA Lewis Research Center 
and is described in detail in References 2.1.16 and 2.1.17. An early version (circa 1983) of this 
code was made available to GE at the start of this project; and it, together with the code 
developed by GE-CRD was installed on a GEAE (GE Aircraft Engines) minicomputer. As 
illustrated in Figures 2 and 3, the computational grids employed by the two codes are radically 
different - it was felt that the results of exercising the two codes for the same propfan case 
would be of interest. 
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Figure 2. Computational Grid • NASPROP-E. 
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Figure 3. Computational Grid - GE-CRD. 



2. 1.2.3 Blade-Loading Parameters 

Blade chordwise loading is input as a normalized function ft, (X), (Equation 2.1.4). This 
is defined as: 


<V 



(X) e 


ik x X 


dX 


Also, 





(X) dX = 1.0 


( 2 . 1 . 7 ) 


Now; from the three-dimensional aerodynamic flow solver we can obtain a chordwise distribu- 
tion of pressure differential Ap = pps - pss as a function of x. 

We thus have: 



Ap (x) dx 


( 2 . 1 . 8 ) 


Write the lift coefficient as 


C, = 


L' 


jPoV r 2 c 


where: 

V r » ao Mr 
c * blade chord. 
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Then 


L' = 


1 

2 r o 


? p, » 2 c c - 


-/ 


te 


Ap (x) dx 


le 


or 



Ap (X) 



— dX 

C L 


Hence, 


f L (x) = Ap(x) - 

7P, v ,c L 


The power absorbed by the rotor can be evaluated from the above as follows: 
At each radial section. 


•' * 2 V ' «L 




( 2 . 1 . 9 ) 


( 2 . 1 . 10 ) 
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The lift direction relative to the direction of rotation is given by 0L (per sketch below). 


L' 



The torque per unit span is thus: 

T = Br L' cos0l 

where, 

B = number of blades 
r = radius of section of interest. 
The power is then given by: 



hub 


Br L' 


cos0^ dr 


tip 

* a J Br "r 2 2 c C l 

hub 


cos0^ dr 
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this becomes 


If we write cos0l = 2 = “ 


Mx __ 1 

tip 


W * 8 * p o 


tiP 

*v/ ( "r “* C L fc) 2 dZ 

hub X 


If it is assumed that Mx is a constant, and wnting 


QD_ 

“t'-zT 

o 


this becomes: 


.tip 

W *?».*.V H «“r */ "r C L fc 2 d2 

hub 


( 2 . 1 . 11 ) 


From the above, it can be seen that the results of the three-dimensional Euler computer 
code can be employed in both the spanwise and chordwise directions, and this is the approach 
followed in the code developed under this contract. However, it is also possible that, bearing 
in mind the long term objective of the development of simplified analytical models for the 
chordwise loading parameters, a hybrid approach could be adopted in which the normalized 
chordwise loading distribution would be obtained from the aerodynamic flow solution as 
described above, while the loading distribution in the spanwise direction would be input 
separately. 

Turbomachinery designers tend to describe blade loading in terms of “swirl”; the 
component of tangential velocity in the absolute frame of reference. It is a measure of the 
work put into the fluid. Propfan references, on the other hand, use “elemental power coeffi- 
cient” in this context These two are related as follows. 


! 
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Power Coefficient: 


c 

p 


p 


P o ^ V 


where 

P = power 
n = rps 
Dt = diameter 
p = density 

Consider an element of fluid passing through 
dP = d(nQ) = nA(rC u )dm 

where 

Cu = Swirl 
and 

dm = poV 0 (2irrdr) 

Assume that, for a propeller, (rCu) in = 0. 
Then, 

dP 2 

= p 0 V 0 2jtQr 2 Q, 


The elemental power coefficient. 


dC P 

dz 


Po v o 2nOr 2 Q, 

Po n °r 
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rotor. The work added to the element is: 



Writing 


U T = 7tnD T , and 0 = — = , 


it can be shown that this equation reduces to: 

® 7l 4 Z 3 M x dz 


2 

7t 3 Z 3 J 



dz 


where J is the advance ratio. 

© is referred to as a swirl coefficient; or, alternatively, a work coefficient. It is input to the 
noise prediction code and used to calculate section lift coefficients as follows. 

Given 



( 2 . 1 . 15 ) 


and 


QA (rCu) = A (aozM-fCu) ; 


writing 
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and assuming zero inlet swirl as before, gives: 


OA (rCJ = *VMr 2 e ; 


and 



Now: 


2n rp V Q 
*o x 

Qr B cosfl^ 


from the above, hence 



2nr p V 

r o X 

c Qr B cosd. 



which yields, after some rearrangement, 

l* 9 . 

cp a 

r o 


where “<x,” the solidity, is given by: 


a = 


cB 

2nr * 
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Thus, 



2z Up 9 
a M r 


(2.1.16) 


(2.1.17) 


The lift coefficients used in the noise code in the evaluation of Equation 2. 1.2 are obtained 
as described above. 

2.1 .2.4 Near Reid Acoustic Model 

In the theoretical development leading to the formulation of Equation 2.1.2, the standard 
far field approximation for R (the source-to-observer distance) was employed. This has the 
effect of lo cating the sound source, for radiation purposes, at the center of the propeller disk 
and is eminently reasonable when the observer is in the acoustic far field. 

In the acoustic near field, however, this compact source assumption is less reasonable; 
individual points on the blades can be considered as acoustic radiators in their own right. The 
inclusion of exact source-to-distance terms in the theoretical derivation greatly increases the 
complexity of the analysis, as can be seen in the fuselage diffraction/refraction formulation of 
Section 23 hereof. 

At the outset of this program, it was felt that the “stationary phase” approach of Hawkings 
and Lowson (Reference 2.1.2) could be employed in this regime. This result is an 
approximation for large values of mB and is valid irrespective of the location of the observer. 
It is thus applicable in the acoustic near field. The result shows that the main contribution to 
the integrals over the blade surface comes from the vicinity of a so-called stationary point. In 
this instance, this is given by: 

R o a o A 3n 

~ 9 r 
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where r s is the stationary radial location, e is the azimuthal location of the source point, and 
Ro is the distance from the center of rotation to the observer. 

In addition, as an alternative approach, it was proposed that, since the aerodynamic Euler 
code calculations essentially evaluate the pressure distribution on a cylindrical surface 
surrounding the rotor, it should be possible to estimate the near field noise directly from the 
flow solution. The main disadvantage to this approach is the necessity of obtaining a solution 
with a sufficiently fine resolution in terms of grid density in the blade-to-blade direction. It 
was felt, however, that with advances in supercomputers, obtaining these solutions in a 
reasonable amount of computing time was not totally infeasible. 

As a preliminary exercise, data/theory comparisons were made using the far field radia- 
tion model already in the computer code. Results were encouraging, leading to the belief that 
development of a sophisticated near field model was, perhaps, not as imperative a task as had 
at first been envisaged. 

It was, thus, decided to adopt the semiempirical approach used by Sulc, et al. (Reference 
2.1.18) for conventional propellers; and to apply such modifications as seemed appropriate 
with regard to current experience. 

Sulc, et al. found that the near field noise of a conventional propeller can be predicted 
reasonably well utilizing a far field formulation (such as, Gutin’s equation) if the source-to- 
observer distance is replaced by the distance between the observer and an equivalent near 
field source. 

In the far field, the range is given by the distance from the observer to the center of rota- 
tion of the propeller, Ro. The equivalent near field range is given in Reference 2.1.18 as: 

Rc = Ro - 0.7 r T 


where 


r T = tip radius. 

Effectively, this is assuming that the propeller noise sources lie on the surface of a sphere, 
centered at the center of the propeller, and having a radius of 0.7 r T . 

However, it may be speculated that the effective source radius will be some function of 
the source frequency (or wavelength) with the effective radius becoming smaller as the 
frequency increases. 

With this in mind, let us write: 

(^) r T 
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where 


x = wavelength of the sound source (propeller); 
K = a constant whose value is to be determined. 


Now: 


2n c 



where 


Mt = tip Mach number. 


If we follow Reference 2.1.18 and state that for a conventional four-bladed propeller with 
a tip Mach number, Mt, of 0.8, 



we have 


0.7 mBHp Q.7 x 4 x Q.8 
S 2n 27t 



Further, assumin g that this “constant” is general, 
„ _ „ 0.7 n 

R = R " nu ^*nr* 

e o mBrU i 
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This shows that as the frequency (represented by the number of blades and the harmonic 
number m) increases, the effective source radius becomes smaller and the effective range 
approaches the far field limit 

The effects of forward flight can be included by adopting the procedure described in 
Reference 2 . 12 . The frequency is multiplied by a factor 

l 


where 


C v = 1 - Mo cos e e 


and 


Mo is the flight Mach number 
e e is the emission angle. 


In addition, the source-to-observer distance, R e , is reduced by a factor C v , resulting in 



( 2 . 1 . 19 ) 


On simplification, this shows that the “near field effect” can be expressed as: 


AdB - S - 
□f 



( 2 . 1 . 20 ) 


For an 8-bladed rotor operating at a tip Mach number equal to the flight Mach number of 
0.8, the near field corrections (calculated at a tip clearance of 0.8 times the diameter), range 
from 1.23 dB at the first BPF (Blade Passing Frequency) harmonic to 0.23 dB at the fifth 
harmonic. 

It is felt that this approach is reasonable but that further study is required to determine 
the “constant” K with a higher degree of confidence. 
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2.1.3 Data/Theory Comparisons (Steady Loading and Thickness Model) 

2.1 .3.1 Aerodynamic 

In the absence of detailed measurements of the aerodynamic flow field around a high 
speed propfan, direct comparisons between the flows calculated by the three-dimensional 
Euler codes and that achieved in real life are not feasible. An indirect comparison, showing 
the spanwise variation of the elemental power coefficient, can be seen in Figure 4 which com- 
pares data from Reference 2.1.19 with values obtained from the results of 3D Euler calcula- 
tions of the flow field of the SR-3 propfan at its design point Before hasty conclusions are 
drawn from the comparison, it should be emphasized that at the time when these calculations 
were made, both flow codes were at a relatively early stage in their development Certain 
anomalies which may be apparent in results presented here have since been corrected. 

Examination of Figure 4 shows the two aerodynamic codes to have considerably higher 
values of elemental power coefficient (especially toward the tip of the blade), than do the 
measured data. An explanation for this is provided in Figure 5, wherein the same spanwise 
information is plotted in the form of swirl coefficient Both aerodynamic codes show finite 
values of swirl at a z-value of unity, implying that the tip of the blade is carrying some load. 
The conversion from swirl coefficient to elemental power coefficient is (from Section 2. 1.2.3, 
Equation 2.1.14) 


9 = (s 3 ? 1 


,dC 
\ — E 
Jl dz 


or, by rearrangement 


_ n 3 z 3 J 
dz 2 


0 . 


From this it can be seen that small differences in swirl coefficient toward the tip of the 
blade will produce (through the Factor z 3 ) much larger discrepancies in dCp/dz than would be 
observed for the same swirl coefficient differential located toward the hub region. 

It was noted during this investigation that calculations of absorbed power by means of the 
equations sketched out in Section 2. 1.2.3 using the swirl coefficients obtained from the 
aerodynamic flow codes did not necessarily produce results in agreement with measurements. 
It was decided to assume that the predicted distribution represented the physics with regard 
to shape, but the levels should be adjusted by a constant factor of SHP mca sured/SHPpredicted 
for each case under consideration. The values of lift coefficient used to calculate the steady- 
loading noise harmonics reflect this adjustment. 
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Figure 4. Comparison of dCp/dz. 
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Figure 5. Propfan Exit Swirl Coefficient Distribution Comparison. 


26 




The second input required from the aerodynamic solution is the normalized chordwise 
load function fL(X), defined in Equation 2.1.10 as: 


f L (X) = 


Ap(X) 

\ p,v?c l ' 


Figures 6 through 8 show examples of this function, which were calculated for the design 
point of the SR-3 propfan by the two 3D Euler solvers available. The presence of two 
NASPROP-E solutions reflects the utilization of a better estimate of blade geometry in the 
solution used to generate Figure 7. It will be noted that, whereas, all three solutions give 
chordwise distributions in the tip region that can be described as rectangular; toward the hub 
there is a divergence from this pattern. This divergence is most pronounced in Figure 8, 
generated from the GE in-house solution. Comparisons with later figures indicate that the 
location and magnitude of the peaks in the loading distribution correspond to the location and 
strength of the shock pattern computed on the blades by the two computer codes. In terms of 
the blading acoustics, it would be anticipated that the rectangular distribution at the tip would 
have greater significance than the distorted distribution at the hub. 

Figures 9 through 12 show contours of relative Mach number on the surface of the SR-3 
blade at its design point, as calculated by the two codes. These, together with Figures 13 and 
14 (which demonstrate the contours in the hub blade-to-blade passage), go some way toward 
explaining the differences in Figures 7 and 8. For example, the supersonic relative flow calcu- 
lated by NASPROP-E in the hub passage is confined to two isolated “bubbles,” whereas in 
the GE calculation, there is a region of supersonic flow extending across the blade passage. It 
is the presence of these large regions of supersonic relative flow in the GE prediction that 
gives rise to the distorted character of the normalized chordwise load distribution. 

One further question concerning the available flow solutions surfaced during this 
investigation. Because of the large amounts of computer time required to obtain a converged 
solution from scratch, the standard operating procedure that was employed was to obtain a 
converged solution at the design point of the propfan, and use this solution as the starting point 
for the various “off-design” cases required. Figure 15 portrays the radial variation of the 
circumferential average of static pressure along the grid lines corresponding to the leading and 
trailing edges of the blading, and extending out into the region beyond the tip of the blades. 
This solution is from the NASPROP-E program calculation of the SR-3 at its design point and 
is regarded as a converged result It will be noted that the limiting value is unity; this is a 
function of the normalization system used in the code. The tip of the blades is denoted by the 
R = 0.5 grid line. Figures 16 through 18 show equivalent plots for different values of helical 
tip Mach number, all of which took the solution of Figure 15 as the starting point. These solu- 
tions were all obtained for constant advance ratio J = 3.06, thus the change in helical tip Mach 
number implies changes in axial and rotational Mach numbers to maintain the value of J. It 
can be seen that, in the region outside the blades, the static pressure “gradient” from blade tip 
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Figure 6. Chordwise Loading - NASPROP-E (1). 
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Figure 7. Chordwise Loading - NASPROP-E (2) 
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Figure 9. Contours of M re i on the Blade Suction Surface • NASPROP-E. 
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Figure 15. Radial Variation of the Circumferential Average of Static Pressure at Blade 
Leading and Trailing Edges - Mn r = 1.15, Moo = 0.8, No Starting Solution. 















39 












40 





to the boundary becomes more and more pronounced as the helical tip Mach number is 
reduced. In contrast, the solution used to produce Figure 19 was run from scratch. Here, the 
region outside the blades appears as uniform as that of Figure 15. 

As a result of these figures, some doubts surfaced as to the validity of the utilization of 
these solutions in the prediction of the noise generated by high speed turboprops. However, 
it was determined that, since the region of acoustic interest was confined to the vicinity of the 
propeller disk, and since all the quantities required were normalized by other parameters from 
the same solution, the use of these results was justified. 

It should be emphasized here that it is not the purpose of this report to criticize or to 
express dissatisfaction with the flow solutions available at the time this work was performed. 
These examples are provided as an illustration of the uncertainties involved in the acoustic 
predictions used in the data/theory comparisons which follow. 

2.1 .3.2 Acoustic 

Comparisons between measured and predicted tone results have been made for three 
model propfans identified as SR-2, SR-3, and SR-6, respectively. These are illustrated in 
Figure 20, which is taken from Reference 2.120. Further theoretical predictions have been 
made for the SR-7L propfan as it performed on the Gulfstream airplane. These are described 
in Section 3.0 of this report As shown in Figure 20, SR-2 and SR-3 each have 8 blades; SR-6 
has 10 blades. SR-2 resembles a conventional propeller in having an unswept planform; 
SR-3 has high sweep and SR-6 falls somewhere between the other two. 

Acoustic and performance data for these propfans were available from References 2.1.20, 
2.121, 2.123 and 2.124. All three model propfans were tested in the NASA Lewis Research 
Center 8- by 6-foot transonic wind tunnel; in addition, data measured from the flights of the 
SR-2 and SR-3 models on the Jetstar airplane were available from Reference 2.1.21. The 
available data points are presented in Table 1. 

Data/theoiy comparisons are presented in different formats in Figures 21 through 31. 
Figures 21 through 24 compare prediction with data for the first six harmonics of BPF at various 
observer angles; while in Figures 25 through 27, the variation of peak BPF tone SPL (sound 
pressure level) with helical tip Mach number at approximately constant C p (power coefficient) 
and J (advance ratio) are shown. These comparisons are discussed in detail below, together 
with possible causes of any discrepancy between predicted and measured levels. 

Figure 21 compares measured and predicted levels for the SR-2 model propfan mounted 
on the Jetstar aircraft. These data (from Reference 2.1.21) were obtained with both boom- 
mounted and fuselage-mounted transducers. Both sets of data (boom and fuselage) were 
taken at the same sideline distance. However, as a result of the different transducer installation 
geometries (namely, the difference in radius of the boom and fuselage), the levels measured 
at equivalent angular locations on the boom and fuselage differed. Consequently, the decision 
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Figure 19. Radial Variation of the Circumferential Average of Sialic Pressure at 
Blade Leading and Trailing Edges - M|jp ~ 0.72, M oo = 0i, No 
Starting Solution. 
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Figure 20. Design Characteristics and Planforms of Various High Speed Propeller 
Models. 
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was made to make the comparisons under free field conditions; and thus, following Reference 
2 i 22, 6 dB were subtracted from the fuselage levels and 4 dB from the boom levels to simu- 
late the free field theoretical results. 


Table 1. Data Used tor Comparison with Theory. 


Propfan 

Reference 

M Ht 

M 0 

a 

J 

SR-3 

2.1.23 

121 

0.85 

1.5 

3.07 



1.14 

0.80 

1.78 

3.07 



1.07 

0.75 

1.79 

3.05 



1.00 

0.70 

1.84 

3.06 



0.863 

0.60 

1.89 

3.05 

SR-3 

2.1.24 

121 

0.85 

1.53 

3.05 



1.14 

0.80 

1.71 

3.06 



1.08 

0.75 

1.79 

3.06 



1.00 

0.70 

1.89 

3.06 



0.93 

0.65 

1.91 

3.06 



0.86 

0.60 

1.91 

3.06 

SR-3 

2.1.21 

1.124 

0.787 

1.998 

3.083 

SR-2 

2.123 

121 

0.85 

1.89 

3.07 

SR-2 

2.121 

1.105 

0.787 

2.094 

3.186 

SR-6 

2.120 

1222 

0.85 

1.02 

3.06 


Two sets of theoretical predictions are compared in Figure 21. This is a result of selections . 
made when the 3D Euler computer codes were being exercised. It was decided at that time 
to exercise the code at a constant advance ratio (J), with varying values of helical tip Mach 
number. In addition, the solutions were obtained for a single value of 33/4, the blade pitch 
angle, for each blade design. This enabled the different flow solutions for each blade to be 
obtained from a single computational grid. For the SR-2 propfan case (Figure 21), flow solu- 
tions were obtained using the NASPROP-E code for a constant advance ratio of 3.06, with 
helical tip Mach numbers of 0.86, 0.95, 1.00, 1.07, 1.15, and 1.20; and the blade pitch angle, 
33/4 of 58°. 

The particular data set plotted in Figure 21 came from a run where the advance ratio was 
3.186, and the helical Mach number (Mht) was 1.105. It was decided to predict this case twice 
using results from the Mht = 1.15 solution (Theory 1), and from the Mht = 1.07 solution 
(Theory 2). Fuselage data were available for six angular locations, with data from the boom 
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Figure 24. SR2—SR6 Wind Tunnel Spectra (130*) Under Free Field Conditions 





Variation of Peak BPF Tone SPL with Helical Tip Mach No. 
for the SR-3 Model Propfan in the NASA Lewis Research Center 
8ft x 6ft Wind Tunnel. Note: WT1 = Ref 2.1.23; WT2 = Ref 2. 



Figure 25. Peak BPF Tone SPL Versus Mhj, Data and Theory (SR-3) Under Free 
Field Conditions. 
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Figure 28. Near Field Calculation Grid, NASPROP-E. 
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above the propfan available for four of these. In general, the predictions agree better with the 
boom data than with that from the fuselage. It is considered probable that, owing to the 
presence of the boundazy layer on the fuselage, the assumption of “flat plate” pressure 
doubling (that is, 6 dB) may be in error. 

Agreement between the boom microphone data and both predictions is generally seen as 
very good; in that, in most instances the two predictions are very close. At certain places, 
however, such as at the third harmonic at e * 78.8, the predictions given by the two codes are 
of the order of lOdB apart. Because in each case the thickness distribution is identical, this 
disparity must come from differences in the calculated spanwise and chordwise loading dis- 
tributions and their interaction with the thickness distribution. 

Figure 22 shows a similar plot for the SR-3 model propfan in flight on the Jetstar aircraft. 
For the SR-3, the same matrix of flow solutions was available with, in this instance, a blade 
pitch angle of 613°. The Jetstar flight was made with a 03/4 of 58.9° and Mht of 1.124. For 
this prediction, the flow solution obtained at Mht = 1.15 was used. Again, the measured data 
were adjusted by subtracting 6 dB from the fuselage microphone values, and 4 dB from the 
boom microphones (where available). Once again, the prediction agrees better with the boom 
data than with the fuselage data, with a major difference being in the fundamental of BPF at 
an observer angle of 69.4°. Here, the prediction is approximately 10 dB higher than the actual 
measurement. 

It is surprising that a discrepancy of this nature should be observed in the fundamental 
BPF tone, especially when the other harmonics are in such good agreement. More reasonable 
is the gradual “drift” between data and theory observed at e = 90.5° with increasing harmonic 
number. It has been noted that the thickness component of the noise prediction tends to hold 
the levels up with increasing harmonic number, especially near the plane-of-rotation. 

In summary then. Figures 21 and 22 demonstrate good agreement with the boom 
microphone data; while data from the fuselage-mounted microphones is consistently lower 
than that measured on the boom, and also lower than the predicted values. This has been 
ascribed to the effect of the boundary layer on the fuselage. 

Figures 23 and 24 compare results obtained in the NASA Lewis Research Center tran- 
sonic wind tunnel with predictions for the SR-2, SR-3, and SR-6 model propfans. These data 
are taken from References 2.120, 2.1.23, and 2.1.24. The prediction program requires, as 
input, local ambient conditions of temperature and pressure (in a wind tunnel environment, 
this amounts to static temperature and pressure), and SHP. For these predictions, tunnel static 
pressure was obtained from Reference 2.1.25, and the temperature and shp were deduced 
from the information in References 2.1.20, 2.1.23, and 2.1.24. 

No flow solutions were available for the SR-6 propfan, owing to problems with the 
NASPROP-E code. Approximations to the spanwise and chordwise loading distributions for 
this propfan were made using data from Reference 2.1.26 for the spanwise variation and a 
simple rectangle for the chordwise distribution. It is recognized that these approximations are 
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less than ideal and will probably account for at least some discrepancies between measured 
and predicted results. The data used for these comparisons were obtained with pressure 
transducers mounted in bleed holes in the ceiling of a perforated-wall wind tunnel working 
section. Where necessary, they have been adjusted to account for a 6 dB error discovered at 
a later date and reported in Reference 2.1.27. Again, the comparisons are made under free 
held conditions. 

Figures 23 and 24 compare measured and predicted spectra at observer angles of 1 10° and 
130°, respectively, at a wind tunnel Mach number of 0.85, giving a helical tip Mach number of 
1.21, at different power coefficients for the three propfan models. The agreement for SR-2 
and SR-3, while not as good as that obtained with the boom microphones on the Jetstar, is 
reasonable. Possibly for the reasons stated above, it is less good for the SR-6 comparison. 

A frequently employed form of presentation (for example, References 2.1.20 through 
2.1.24) is to plot the variation in peak Blade Passing Frequency SPL as a function of the blade 
helical tip Mach number, Mht. at constant advance ratio, J. Figures 25 through 27 
demonstrate this variation for the SR-3 propfan as tested in the NASA Lewis Research Center 
8- by 6-foot transonic wind tunnel and on the Jetstar airplane. In all cases, the theoretical 
predictions presented in these figures were made under the measurement conditions (but free 
field), and such adjustments as were made (with the exception of conversion to free field) were 
applied to both sets of information. 

A major difference between the data gathered in the wind tunnel and that obtained on the 
Jetstar is found in the ambient pressure at which the measurements were made. Jetstar data 
were primarily taken at 30,000-feet altitude, with constant ambient pressure. In the wind 
tunnel, a change in tunnel Mach number is accompanied by a corresponding change in static 
(acoustic ambient) pressure; consequently, the tone levels shown in Figure 25 have been 
adjusted using 

SPL(plot) = SPL - 20*logio(Pamb/Pref) 
where P re f = 11.11 psi, the tunnel static pressure at M 0 = 0.8. 

Likewise, in Figure 26, with the exception of the points at Mht = 0.751, all the data were 
taken at constant altitude (30,000 feet) and, hence, constant pressure. The values measured 
at Mht — 0.751 were adjusted to reflect the change from 20,000 to 30,000 feet. Additionally, 
the boom data were adjusted to free field by subtracting 4 dB, as before; while 6 dB were 
subtracted from the fuselage readings. Finally, Figure 27 combines Jetstar and wind tunnel 
data and predictions (for this figure, the wind tunnel conditions of Figure 25 were taken as 
standard). Two adjustments were applied to the Jetstar data of Figure 26; namely: 

1. From 4364 — > 11.111 psi .... + 8.12 dB (ambient pressure) 

2. From 2.65 — > 4.01 feet ..... - 3.60 dB (measurement distance) 

Total adjustment = + 4.52 dB. 
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Examina tion of these figures shows that the measured and predicted values of Figure 25 
are in good agreement, up to a tip helical Mach number of about 1.1. Beyond this point, the 
predicted levels continue to increase with increasing tip helical Mach number, but the 
measured data appear to level off. This flattening of the curve is not seen in the Jetstar data 
(Figures 26 and 27), wherein the boom data can be seen to be in good agreement with the 
(wind tunnel) predictions, but the fuselage microphones show a consistently lower level. 

In Reference 2.1.28, Dittmar reexamined the SR-2 and SR-3 data from the wind tunnel 
and Jetstar. He concluded that whereas both exhibited a tendency for the peak SPL in the 
blade passing frequency tone to remain constant at high helical tip Mach numbers; the flat- 
tening out occurred earlier with the SR-2 than with the SR-3. The SR-3 propfan was used for 
the present study because of the availability of a greater number of flow solutions. However, 
the agreement between data and theory as illustrated in Figures 21, 23, and 24 for the SR-2 
propfan is such as to suggest that the theoretical prediction would exhibit the same trend as 
the data for this unswept planform. 

In conclusion, these results show that the approach of combining nonlinear aerodynamics 
with linear acoustics has led to good agreement with data acquired with the boom microphones 
on the Jetstar for both the SR-2 and SR-3 model propfans. Of the three measurement “envi- 
ronments’* available, the boom microphones most closely model the free field assumed in the 
theory. Data from the Jetstar fuselage is thought to be affected by the fuselage boundary layer 
(Reference 2.122), and the different trends with increasing helical tip Mach number observed 
on the Jetstar and in the wind tunnel are confusing. The SR-6 results shown suggest that the 
assumptions made as to spanwise and chordwise loading distributions, although reasonable, 
were not as good a representation of the flow field on the blading as that obtained from the 
Euler solutions for SR-3 and SR-2. The agreement for the SR-3 and SR-2 propfans suggests 
that the worries expressed over the degree of convergence of the flow solution outside the 
blading were unfounded, and also that flow solutions obtained at operating conditions that are 
not identical to those under which acoustic data were taken can still provide useful and valid 
input to an acoustic calculation. 

The predicted results used in the above-described data/theoiy comparisons were obtained 
from a far field, frequency-domain, acoustic model with a semiempirical adjustment to take 
near field effects into consideration. A second approach to this problem is discussed briefly 
in Reference 2.128 and in Section 2.1.12 of this report, and is given here in more detail. The 
basic premise is that, since the three-dimensional flow solution describes the blade-to-blade 
pressure field in a frame-of-reference (locked to the rotor), by allowing this field to rotate, the 
time-varying pressure perceived by an observer is obtained. This can then be Fourier analyzed 
to obtain the acoustic near field. 

A NASPROP-E solution of the design point case (Mo = 0.8, Mht * 1.15) of the SR-3 
propfan with 24,750 grid points (21 in the blade-to-blade direction) was utilized to investigate 
the feasibility of the use of 3D Euler flow solutions in the prediction of the acoustic near field 
of a high speed propfan. 
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A projection of the first meridional surface was used to establish the radial and axial loca- 
tions at which the calculation should be performed. These were selected to coincide with the 
angular positions of the microphones in the NASA 8- by 6-foot wind tunnel test (Reference 
2.1.23). The selected positions are identified in Figure 28. 

For the current investigation, a Fourier analysis of the rotating pressure field was con- 
ducted at two values of (R/Rtip) for each angular position under consideration (75°, 90°, 100°, 
and 110°); namely R/Rtip = 1.0, in all cases, and R/Rtip = 2.7965, except at 75°, where the 
second calculation was performed at R/Rdp = 2.18. At each location, the computed pressure 
distribution in the blade-to-blade direction was interpolated to 101 points and Fourier 
analyzed to extract the harmonic content Adjustments were made to these results to account 
for differences in static pressure between the wind tunnel and computational conditions, and 
differences between the computed and measured power coefficients under which the propfan 
was operating. The computed results were then extrapolated out to the sideline distance at 
which the measurements were made, using a simple 1/r 2 law. 

Figures 29 through 32 compare these data with the results obtained as described above. 
The data used for comparison were obtained from Reference 2.1.23 and have been adjusted 
to free field conditions (-6 dB), and corrected for an error ( + 6 dB) discovered after the original 
report was published (Reference 2 . 127 ). These corrected results are best examined in light 
of Figure 28, which gives some indication of how the computational grid used may influence 
the calculation. 

At 75°, there is only one data point for comparison (Figure 29). The fundamentals of BPF 
calculated at the two radii are in close agreement, but are approximately 10 dB below the 
measured value. Higher harmonics fall off more rapidly at the greater radius. Consideration 
of Figure 28 shows that at R/Rtip = 1, the 75° ray is in a region where the axial grid is in tran- 
sition from the coarse spacing used over the leading edge of the nacelle, to the much finer 
spacing used on the blade. At R/Rtip = 2.18, the grid is considerably coarser. It appears that 
the pressure field in both locations is affected by the grid, with the greater refinement of the 
points nearest the blade contributing to the slower falloff with increasing harmonic number. 

Similar arguments can be applied at 90° (Figure 30). An apparently very good agreement 
of the fundamental BPF measured with that calculated at the tip radius is evident. The falloff 
of the calculated harmonics suggests that greater resolution is required in the blade-to-blade 
direction. 

At 100°, the story appears somewhat different (Figure 31). Here the values calculated 
using the nearer radius are actually higher than the measured data, while the fundamental BPF 
calculated from the further radius is in very good agreement with measurement. This is 
believed to be due to the fact that at R/Rtip = 1.0, the calculation is performed at the physi- 
cal tip of the blade. It is believed that there is a genuine near-field effect in this calculation, 
whose rapid decay rate is not accounted for by the simple distance rule applied to the calculated 
values. 
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At 110° (Figure 32), the grid at the outer radius is extremely coarse. This is believed to 
be the reason for the very low values calculated. At the inner radius, on the other hand, agree- 
ment with data is obtained for the first three harmonics of BPF. This is also the angle at which 
the highest levels of noise were measured, so agreement between data and theory here is 
encouraging. Again, the degree of agreement in the higher harmonics is thought to be grid- 
dependent 

In view of the difficulties involved in interpreting these predictions without the benefit of 
measured data, and bearing in mind the added computational time and cost of obtaining 
solutions with a sufficient number of grid points in the blade-to-blade direction to guarantee 
adequate resolution for the acoustic calculation, this approach was abandoned until such time 
as available computers are found to have decreased in cost and increased in speed; each by an 
order of magnitude. 

2.1.4 Conclusions 

The following conclusions were drawn: 

• A far field, frequency domain, noncompact source, linear acoustic model has been 
developed and programmed to predict the steady loading and thickness com- 
ponents of the noise generated by a high speed single rotation propfan. 

• This model has been combined with loading distributions obtained from three- 
dimensional nonlinear aerodynamic Euler code flow solutions to predict the noise 
generated by two different model propfans, both in flight and in the wind tunnel; 
simplified loading distributions have been used to predict the noise from a third 
propfan. 

• Comparisons between measured and predicted data show, on the whole, good 
agreement; differences between data and theory are of the same order as differ- 
ences between data sets. 

• The results show that nonconvergence of the flow solution in the region outside 
the blading does not necessarily compromise the acoustic calculation. 

• The use of flow solutions directly in the calculation of near-field noise requires 
greater grid resolution than is practical at this point in time. 

2.2 Installation Effects on Single Rotation Propeller Noise 

2.2.1 Introduction 

The objective of this study was to develop an approximate, analytical model to compute 
the additional noise produced by a single rotation propeller as a result of installation effects. 
The three installation effects considered herein are: the angle-of-attack (defined as the angle 
between the freestream and the propeller axis), the effect of wing-induced flow, and the effect 
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of flow induced by an axisymmetric body of revolution representing the fuselage. The calcu- 
lation procedure is as follows. 

First, the steady loading and thickness noise of the propeller is calculated. Next, the non- 
uniform axial and tangential flow fields induced at the propeller plane are developed for the 
three installation effects to be evaluated. Quasi-steady propeller aerodynamic theory is used 
to compute the induced unsteady axial and tangential forces on the propeller due to these non- 
uniform flow fields. In order to better account for unsteady effects, a phase lag deduced from 
linearized, two-dimensional, compressible, flat plate gust response theory is used in conjunc- 
tion with the quasi-steady analysis to estimate the unsteady forces, which are then utilized in 
conjunction with an acoustic analysis to derive the far field noise. 

The noise field due to installation effects has been found experimentally to be distinctly 
nonaxially symmetric (with respect to the propeller axis), and hence, the phase between the 
noise field due to steady loading (and thickness) and that due to unsteady loading must be 
accounted for. 

The following material is discussed in four categories: 

• Calculation of nonuniform flow fields due to angle-of-attack, wing loading, and an 
axisymmetric body of revolution representing the fuselage 

• Unsteady force calculation 

• Acoustic calculation linking the far field noise due to unsteady loading, steady load- 
ing, and thickness 

• Data/theory comparisons. 

2.2.2 Calculation of Nonuniforni Row Reids 
2 .2.2.1 Angle-of-Attack 

This source of nonuniformity is obviously explicitly calculable and yields the simple result 
that, with the coordinate system shown in Figure 33, for small angles-of-attack (a) such that 
cosa - 1, the angle-of-attack induces a tangential velocity perturbation (v): 

v = UsinaCOs<j> 

where U, a, and $ are defined in Figure 33. 

2J22.J2 Wing Loading 

The flow field due to wing loading is calculated according to a model discussed in Article 
12.4 of Reference 2.2.1. In this model, the flow field is associated with an airfoil of span 2s' 
(with 2s denoting the true wing span) with constant circulation K and the associated trailing 
vortex system (the reader is referred to Figure 93 of Reference 2.2.1). 
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Figure 33. Definition Sketch of Angle-of-Attack and r-$ Coordinate System. 


For an airfoil with elliptic loading, (s'/s ) can be shown to be (ir/4), and, given the lift 
coefficient and aspect ratio of the wing (Cl, AR), we can calculate K by: 

K = C L n 3 AR (2.2.2) 

Knowing K and s' , the full expression for the three-dimensional flow field due to the above 
horseshoe vortex system is given on Pages 158 and 159 of Reference 2.2.1. These expressions 
are used to calculate the induced axial and tangential velocities in terms of radial location, r, 
and azimuthal angle, <t>. Apart from Cl, and s, an input to such a calculation (in addition to r, 
and <>) is the location of the center of the propeller disk relative to the wing. 

2.2.2.3 Row Reid Due to an Axisymmetric Body of Revolution 

The axisymmetric body of revolution chosen is a Rankine solid (Reference 2.2.2). This is 
the closed streamline shape generated by superposition of a uniform flow (U), a point source 
and sink (of strength “m”) separated by a distance “2a” in the direction of the stream. If the 
total length of the body is “21,” and its maximum radius “h”, then m and a are determined by 
the solution of: 

4a lm a U(l 2 -a 2 ) 2 (2.2.3) 


4ma a Uh 2 V h 2 -a 2 (2.2.4) 

Thus, knowing the total length of the body and a slenderness ratio parameter (such as 
“h/1”), its flow field can be computed. Finally, the chosen body of revolution must be located 
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relative to the wing. This can be accomplished, for example, by specifying the location of the 
forward stagnation point of the body relative to the wing. 

2 .2 .2.4 Total Row Reid 

Given the above analytical expressions for the three sources of nonuniform flow, the resul- 
tant flow field is computed at several fixed radii as a function of azimuthal angle, <t>, in the 
propeller disk plane (Figure 34). At each fixed radius, the axial and tangential nonuniform 
velocity components are developed in a Fourier series in $ of the form: 

N 

^ V 01 ( “(♦-♦],)) 

k=l 

”N” is chosen typically as 7 or 8. This Fourier decomposition is carried out numerically using 
an FFT (Fast Fourier Transform) algorithm. 



Figure 34. Circles Around Which Fourier Analysis of u, v Distortions is Carried Out. 


2.2.3 Unsteady Force Calculation 
2.2.3. 1 Quasi-Steady Formulation 

Let u(r, ■»), v(r, $) denote externally contributed nonuniform axial and tangential velocities 
incident on the plane of the propeller disk, “u” is taken positive in the downstream direction 
and “v” is taken positive in the direction of propeller rotation. 

If “R” denotes the tip radius of the propeller disk and “Rh” the hub radius, u and v are 
defined for Rh £ r* R, and 0 s <j> 2 2tt. 
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At any radius, the quasi-steady forces (thrust and tangential force) are obtained in the 
functional form: 

F = n 2 f(x) 

where “n” denotes the propeller rotational velocity in radians per second, and “x” the speed 
ratio “U/nR,” where U is the advance velocity of the propeller. 

The change in the force, F, due to u and v can be obtained as follows. 

In the quasi-steady approximation, changes in n and x due to u and v, will be: 

5 Q = • - 
r 

The resulting change in F can be written: 

5F = 2fif(A)6fi + (§f ) 

* [-♦>(!)] h 

Thus, if u, v are nondimensionalized by introducing: 



The quasi-steady theory is needed to obtain F and aF/ax . 

Now, aF /ax can be approximated by a central difference formula: 

9F F(A+e) - FfX-e) 

3X ~ 2e 

where “e” is a small increment. 

It can be seen that the theory is thus only required to obtain the force F for arbitrary x 
since a F/ax can be approximated by numerical differentiation as outlined above. 
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The theory needed to compute both axial and tangential forces is based on that of Glauert 
as presented in Reference 223 . The formulation is based on incompressible flow theory and 
the inputs needed are: relative radius (r/R) of interest, solidity (chord/transverse spacing) at 
this radius of interest, number of blades, specification of the section lift and drag coefficients 
as functions of angle-of-attack (including the orientation of the reference line with reference 
to which the angle-of-attack is specified), and the speed ratio, x. = U/n R. The version used 
for the work described herein incorporates the Prandtl tip loss factor. (There is a misprint in 
the development in Reference 2 23 - Equations 5.5 and 5.6 of Chapter VTI - accounting for 
the tip loss factor. This misprint is also pointed out by EJE. Larrabee in SAE Paper 790585.) 
The essence of the analysis is the solution of two simultaneous nonlinear algebraic equations 
(for the so-called axial and tangential interference factors.) 

22.3.2 Phase Lag Effects 

The quasi-steady theory discussed above is felt to be the most appropriate for the large- 
scale low-frequency distortions resulting from the installation geometries addressed in Section 
2.2.2. It is accepted that, in a quasi-steady theory, phase lag or lead effects are not accounted 
for. On the other hand, features of the problem related to the fact that a propeller (possibly 
with a large number of blades) is involved are recognized. To improve the accounting for 
unsteady aerodynamic effects, phase lag (or lead) effects have been included in the develop- 
ment, based on a two-dimensional, linearized, compressible flow gust response theory. The 
expressions used are those of Goldstein (Reference 22.4). At any given radius, if the rotor 
angular velocity is n , the relative velocity V r is taken as: 

VU 2 + n 2 r 2 , 

and the relative Mach number, M r , is just the same quantity divided by the ambient speed of 
sound. The reduced frequency, a*, is a function of the Fourier component of the relevant non- 
uniform velocity and can be written: ok = knc/2V r where “c” is the propeller blade chord at 
radius “r.” Defining a parameter, p, as ckMr/(l-M r ) 2 Goldstein’s Equation 3.71 of Reference 
2.2.4 is used if p is less than unity, and Mr < 1 . Otherwise Equation 3.70 of the same reference 
is used. It should be emphasized that, in the present study, only the phase lag (that is, the argu- 
ment of the complex numbers in Equations 3.70 or 3.71 of Reference 22.4) is deduced from 
the linearized, two-dimensional gust response theory. The amplitude employed is computed 
from quasi-steady theoiy. The data/theoiy comparisons of Section 2.2.5 show the effects of 
including phase lag in the calculations. 

2.2.4 Acoustic Calculation 

In the work described herein, the acoustic source distributions are assumed to be acous- 
tically compact in the chordwise direction. In other words, a line source model is employed. 
At time “t = 0,” let the location of this line source be as shown in Figure 35, and let the sweep 
be as shown in Figure 36. Thus, “ai(r)” and “Ax(r)" define the radial lean and axial sweep of 
the blades. The following results were derived for the purposes of this work by the use of axial 
Fourier transforms and the method of stationary phase, but they are in agreement with results 
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Figure 35. Definition Sketch for Radial Lean. 



Figure 36. Definition Sketch for Axial Sweep. 


given by Hawkings and Lowson (Reference 2.2-5) and Hanson (Reference 2 . 2 . 6 ), and hence, 
only the final result of the analysis is shown. 


The contribution to “p 0 ” from the blade element lying between “r” and ”r + dr” is given 
by: 


dp o inflgfl , [ nBflr sin 9 1 

d T = 2nRc (l-MpCoiQ ) ±nB [ c (I-Mpcose) J 

[j [ fg^SSSr 1 <♦-!> * nBo i (r) • (i-S'Jhsr]] 

r f cos0 cf. jnBp cQhl 1 

zo - 9°. + 0 ij 

|_ ( 1-Mj. cose ) + Qr ( 1-Mj.cose)* J 


(2.2.5) 


In Equation 2 2 * 5 , the far field noise is expressed in terms of R, e (the coordinates of the 
observer relative to the “retarded” propeller disk location.) The relationship between these 
coordinates and Rc, e c (the coordinates of the observer relative to the current propeller disk 
location) is shown in Figures 37 and 38. In Equation 2.2.5, wherever two signs appear, the 
upper sign refers to counterclockwise rotor rotation, and the lower sign to clockwise rotation 
(forward looking aft). 

Now consider the noise due to unsteady loading. Let the unsteady loading per blade per 
unit spanwise length in a frame of reference fixed with a rotating blade be expressed in the 
form: 


Unsteady z component 


L f k cos(kflt) 
k=l 


( 2 . 2 . 6 ) 


and 


Unsteady t component * J f 

k=l 


♦k 


cos(kftt) 


(2.2.7) 


Both fzk and fyk are positive functions of “r.” Also, let Pzk(r) and p«j>k(r) denote the 
azimuthal angles in the coordinate system of Figure 33 where these unsteady components 
achieve their maximum positive values. Then, unsteady loading noise from blade elements 
lying between “r” and “r + dr” contributes to the far field pressure a quantity dpk/dr, given by 
the sum of: 
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where 


/ dp k\ inB 2 0 . r nBQrsinQ 1 

^ dr y 4tiRc (l-MpCosS ) * 7±nB+k) ^(l-MpCosO) J 

•*f [•> [ S (1-!£cm 0 ) * (± “ B ’ k) (♦•!) 1 ,lBo l (r) * ( l-M^eosfl ) ] 

exp ( -jkP zk ) f 2k cos0 _ c exp ( -jkp Ak ) f 6k ( nB±k) 

(1-Mp co30) * nBfir 


/ dp k\ inB 2 0 _ 

. nBfirsin0 1 

y dr 4 tiRc (1-Mpcos0) 

' J (±nB-k) c (l~Mj.cos0 ) J 


“» J [*(£“«$ * C1 “ i * k) (♦•§)* DB0 I (r) * ( l-Spos«) 


««P (JfcP zk ) f zk c °* 9 
( 1-Mj.cos0 ) 


+ nB Or 


(qB + k) 


The total far field acoustic pressure can now be computed as: 



(2.2.9) 


( 2 . 2 . 10 ) 


where rh and rt denote the inner and outer radii of the propeller disk, respectively, and where 
“p” is a complex number whose amplitude gives the amplitude of the far field pressure in the 
nth harmonic of blade passing frequency. 

The above formulation could, perhaps, be expressed more concisely with more extensive 
use of complex notation. However, it was felt that it was preferable to restrict consideration 
to positive n and k and, also, to consider both clockwise and counterclockwise rotating 
propellers. These interests, and the desire to be accurate concerning phase relations, have led 
to the above formulations. One min or addition has been made to the above equation. 

The existence of unsteady thickness noise has been recognized for the following reason. 
The thickness noise arises from a source term related to the product of the relative velocity 
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incident on a blade section and its thickness. Since this relative velocity varies in magnitude 
around the propeller disk, unsteady thickness noise can arise. Consequently, this source has 
been taken into account in the computer program associated with this work. 

2.2.5 Comparisons Between Data and Predictions 

Data/theory comparisons have principally been carried out using data from Reference 
22 . 1 . Only angle-of-attack effects have been examined in this case. Predictions have also been 
made for the SR-3 propeller data of Reference 22.8. Assumptions made in order to perform 
the comparisons are noted below: 

• Blade section lift and drag coefficients (versus angle-of-attack) were calculated 
from performance maps given for the SR-2 model propeller in Reference 2.2.9, 
and for the SR-3 model propeller in Reference 2.2.10. Model propeller blade 
geometry information was also obtained from these sources. 

• The SR-2 noise data of Reference 22.7 were obtained at both an 8° and 10° angle- 
of-attack. The theoretical calculations were performed at an angle-of-attack of 9°. 

In addition, the SR-2 tested acoustically and described in Reference 2.2.7 differs 
from that utilized for the performance evaluation of Reference 2.2.9; in that, due 
to power limitations, a four-bladed configuration was used for the acoustic data 
acquisition, rather than the eight-bladed configuration of Reference 2.2.9. It was 
assumed that this change did not afreet the behavior of the lift and drag versus 
angle-of-attack characteristics. 

• The representative airfoil thickness was taken as one-fourth the maximum thick- 
ness at a given radial location. 

• The measurement angles quoted in the acoustic data are taken to observe angles 
as shown in Figure 38. This neglects shear layer corrections in the data of 
Reference 22.7 (which should be small in view of the low axial velocities of the 
freejet). 

• Because the phase lag formulae discussed earlier apply only to convected two- 
dimensional gusts, their use in the present context could be challenged. There- 
fore, results are presented both with and without the use of such phase lag. 

Comparisons for the SR-2 low tip speed case are contained in Figures 39 through 47. These 
comparisons are rather encouraging (especially in the case of Figures 46 and 47), with regard 
to the data scatter. The SR-2 high tip speed comparisons are also encouraging (Figures 48 
through 56), although clearly both the measured and predicted effects of angle-of-attack are 
small and hence the comparisons are somewhat inconclusive. 

Some theoretical calculations relevant to the data reported in Reference 2.2.8 are shown 
in Figures 57 through 59. Figure 59 is interesting in that it shows that, if the phase lag in the 
force response is neglected, the theoretical predictions are in accord with the comments in 
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Figure 39. Four-Bladed SR-2, 443-fps lip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Noise at Zero Angle-of-Attack Versus 0 at 
35-inch Radius from Center of Propeller. 



Figure 40. Four-Bladed SR-2, 443-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increased Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 8. 4> = 0*, Phase Lag in Force Response 
Neglected in Theoretical Prediction. 
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Figure 41. Four-BIaded SR-2, 443-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22 . 7 , Symbols are Measured 
Data, and Full line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle- 
of-Attack (Relative to Zero Angle-of-Attack) Versus 0. <J> = 0°, Phase Lag in Force Response 
Included in Theoretical Prediction. 



Figure 42. Four-BIaded SR-2, 443-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9* Angle- 
of-Attack (Relative to Zero Angie-of-Attack) Versus 0. <J> = 90 s , Phase Lag in Force Response 
Neglected in Theoretical Prediction. 
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Figure 43. Four-Bladed SR-2, 443-fps Up Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22 . 1 , Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angie- 
of-Attack (Relative to Zero Angle-of-Attack) Versus ©.<(> = 90 s , Phase Lag in Force Response 
Included in Theoretical Prediction. 



Figure 44. Four-Bladed SR-2, 443-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22 . 1 , Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/ Decrease of Noise Due to 9 s Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 0. <|> = 270°, Phase Lag in Force Response 
Neglected in Theoretical Prediction. 
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Figure 45. Four-Bladed SR-2, 443-fps Up Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angie-of-Attack) Versus 0. 4> = 270*, Phase Lag in Force Response 
Included in Theoretical Prediction. 



Figure 46. Four-Bladed SR-2, 443-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/ Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus <f>, 0 = 90°, Phase Lag in Force Response 
Neglected in Theoretical Prediction. 
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Figure 47. Four-Bladed SR-2, 443-fps Up Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angie-of-Attack) Versus <J>, 9 = 90*, Phase Lag in Force Response 
Included in Theoretical Prediction. 
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Figure 48. Four-Bladed SR-2, 745-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental 

Blade Passing Frequency Noise. Acoustic Data Source is Reference 22 . 7 , Symbols are 
Measured Data, and Full Line is Theoretical Prediction. Noise at Zero Angie-of-Attack 
Versus 9 at 35-inch Radius from Center of Propeller. 
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Figure 49. Four-Bladed SR-2, 745-fpa Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/ Decrease of Noise Due to 9 s Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 0. <J> * 0°, Phase Lag in Force Response is 
Neglected in Theoretical Prediction. 



Figure 50. Four-Bladed SR-2, 745-fps Tip Speed, 100- fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angie-of-Attack) Versus 0. 4> = 0°, Phase Lag in Force Response is 
Included in Theoretical Prediction. 
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Figure 51. Four-BIaded SR-2, 745-fps Up Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angie-of-Attack) Versus 9. <f> * 90 s , Phase Lag in Force Response is 
Neglected in Theoretical Prediction. 



0C 

Figure 52. Four-Bladed SR-2, 745-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angie-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 0. <j> * 90°, Phase Lag in Force Response is 
Included in Theoretical Prediction. 
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Figure 53. Four-Bladed SR-2, 745-fps Tip Speed, 100- fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 6. <J> = 270°, Phase Lag in Force Response is 
Neglected in Theoretical Prediction. 



Figure 54. Four-Bladed SR-2, 745-fps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9 s Angle-of- 
Attack (Relative to Zero Angle-of-Attack) Versus 9. <J> = 270 s , Phase Lag in Force Response 
is Included in Theoretical Prediction. 
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Figure 55. Four-Bladed SR-2, 745-fps Up Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 22.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9° Angie-of- 
Attacfc (Relative to Zero Angler-Attack) Versus 4>, 0 = 90*, Phase Lag in Force Response 
is Neglected hi Theoretical Prediction. 



Figure 56. Four-Bladed SR-2, 745-fjps Tip Speed, 100-fps Forward Flight Speed, Fundamental Blade 
Passing Frequency Noise. Acoustic Data Source is Reference 2.2.7, Symbols are Measured 
Data, and Full Line is Theoretical Prediction. Increase/Decrease of Noise Due to 9* Angle-of- 
Attack (Relative to Zero Angler Attack) Versus <J>, 8 =* 90*, Phase Lag in Force Response 
is Included in Theoretical Prediction. 
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Figure 57. SR-3 Propeller, Forward Flight Mach Number = 0.8, Advance Ratio J = 3.06. 

Theoretical Predictions for Increase/Decrease of Noise Due to 4° Angle-of-Attack 
(Relative to Zero Angie-of-Attack) Shown by Symbols. Open Symbols are for 4> = 0° 
(North Wall of Reference 2.2.8), and Closed Symbols are for <j> = 180 ° (South Wall of 
2.2.8). Phase Lag in Force Response Neglected in Theoretical Prediction. 
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Figure 58. SR-3 Propeller, Forward Flight Mach Number = 0.8, Advance Ratio J = 3.06. 

Theoretical Predictions for Increase/Decrease of Noise Due to 4 s Angie-of-Attack 
(Relative to Zero Angle-of-Attack) Shown by Symbols. Open Symbols are for 
4 , = 0° (North Wall of Reference 2.2.8), and Closed Symbols are for 4> = 180 s (South 
Wall of Reference 2.2.8). Phase Lag in Force Response Neglected in Theoretical 
Prediction. Sweep of SR-3 Propeller Artifically Suppressed. 
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Figure 59. SR-3 Propeller, Forward Flight Mach Number = 0.8, Advance Ratio J = 3.06. Theoretical 
Predictions for Increase/ Decrease of Noise Due to 4* An gle-of- Attack (Relative to Zero Angie' 
of-Attack) Shown by Symbols. Open Symbols are for <£ = 0° (North Wall of Reference 2.2.8) 
and Closed Symbols are for $ = 180° (South Wall of Reference 2.2.8). Phase Lag in Force 
Response Included in Theoretical Prediction. 
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Reference 2.2.8; namely, that due to angle-of-attack, the north wall data consistently showed 
an increase in noise greater than that seen on the south wall (which sometimes showed a noise 
decrease). 

It should be noted that the results presented in Figures 57 through 59 are plotted against 
“observer” rather than “emitted” angles. (The large forward flight Mach number of 0.8 leads 
to significant differences between current and retarded coordinates). Figure 58 shows that, if 
we artificially neglect the large sweep of the SR-3 propeller, the north-south wall asymmetry 
is predicted to disappear. A comment at the end of Reference 2.2.8 attributes the observed 
asymmetry in noise to such sweep, and the comparison of Figures 57 and 58 supports this view. 
Figure 59 shows that inclusion of a phase lag formula for convected two-dimensional gusts 
apparently gives an asymmetry exactly the reverse of that seen in Figure 57. This also is 
contrary to the data reported in Reference 2.2.8. 

2.2.6 Concluding Remarks 

The data/theory comparisons of Figures 39 through 56 are believed to constitute a some- 
what stringent evaluation in that, essentially, differences in decibels due to angle-of-attack at 
specific microphone locations are being compared. These differences were obtained from 
measured, tabulated data which (in the case of Figures 39 through 56) were kindly supplied 
by the author of Reference 22 . 1 . Less stringent tests would involve some smoothed version 
of the data, perhaps the examination of absolute levels in the form of contour plots, rather 
than specific differences due to angle-of-attack at precise microphone locations. 

As illustrated (Figures 39 through 59), the areas of discrepancy are, by and large, confined 
to angles close to the propeller axis. Agreement between data and prediction improves 
significantly around the plane of rotation. Comparison between propeller noise prediction 
and data at angular locations close to the axis of rotation is always somewhat difficult, due to 
the low levels of noise present at such locations. The data shown in Figures 46 and 55 are of 
special interest, because there is some redundancy in the data shown here, and hence, some 
idea can be gained of the degree of data scatter. The results depicted in Figures 57 and 59 
definitely support a conclusion that a phase lag formula deduced from two-dimensional, 
convected gust response theory for flat plate airfoils should not be used at this time. 

It is recommended that the present analysis, without the above-mentioned phase lag term, 
be used to determine the installation effect on single rotation propeller noise. The result of 
such an exercise should be useful, particularly for angles around the plane of rotation. 

2.3 Scattering of Propeller Noise by Aircraft Fuselage and Wings 

2.3.1 Introduction 

When free-field predictions of propeller noise are compared to measured flight and/or 
wind tunnel data, it is recognized that transducers installed on fuselage surfaces, or on wind 
tunnel walls do not record free-field acoustic radiation levels (see, for example. References 
23.1 and 23.2). Several corrections become necessary in order to account for various forms 
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of acoustic scattering, such as reflection, refraction, and diffraction. Quantitative estimates of 
these effects are necessary, not only for comparison of predicted and measured propeller noise, 
but also for assessing the possible advantages in engine placement relative to fuselage and wing 
surfaces. Further, the nature of the sound field and its spatial distribution on and around the 
fuselage is of great relevance in the design of fuselage structures for minimizing sound trans- 
mission into the cabin. 

The next subsection presents the nature of the problem, with respect to reflection, refrac- 
tion, and diffraction. This is followed (in Section 2 32 ) by a detailed discussion of the analysis 
of propeller noise scattering by a cylindrical fuselage. The incident field is represented by pairs 
of helical waves emitted upstream and downstream of the propeller disk. The frequencies of 
the emitted field are determined by the harmonics of propeller blade passing frequency, and 
the associated spinning mode numbers are equal to the number of propeller blades multiplied 
by the relevant harmonic number. The radial dependence of the sound field is represented 
by cylindrical Hankel functions of order equal to the spinning mode number. Scattering calcu- 
lations are made for the fundamental and two higher harmonics of the blade passing frequency. 
Results are presented to show the increase and decrease of the sound field on and around the 
fuselage surface relative to the incident field as a function of relevant parameters. 

2.3.1. 1 Nature of the Problem with Respect to Reflection, Refraction, 
and Diffraction 

23.1.1.1 Reflection 

a. Rigid Flat Surface • When a plane wave in an ambient homogeneous medium is 
incident upon a flat rigid surface of infinite extent, it will be reflected with a reflec- 
tion coefficient of unity and in phase with the incident field. The combined field 
results in a pressure doubling; that is, a 6 dB increase relative to the incident field. 
This pressure doubling is independent of the angle of incidence, except at grazing 
incidence. 

b. Nonrigid Flat Surface - Under the same conditions as a rigid flat surface (except 
that the surface is now being replaced by a locally reacting, nonrigid surface), the 
amplitude and phase of the reflection coefficient now become dependent, not only 
on the surface admittance, but also on the angle of incidence. The amplitude of 
the combined field (incident plus reflected) is no longer double that of the inci- 
dent field; it is lower. The reduction relative to 6 dB depends on the amplitude 
and phase of the reflected field. 

c. Rigid Curved Surface • Consider a condition similar to (a) above (except that the 
surface is no longer flat, but is rather curved in the form of a cylindrical fuselage). 
Acoustic focusing or defocusing may now be expected. If the surface is curved away 
from the side where the incident sound is coming from (that is, convex), defocusing 
produces a surface sound pressure increase of less than 6dB. Such a defocusing 
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effect was measured by Hubbard (Reference 2 .3.3) and is recommended by the 
SAE for estimating reflection from a circular wall. 

In addition to the defocusing effect, this configuration will also have an acoustic 
shadow region away from the side where the sound field is incident. In fact, sound 
can reach the shadow region as a result of diffraction. The actual level on the 
fuselage surface in the shadow region depends on the frequency, being lower for 
higher nondimensional frequencies based on the fuselage radius. 

d. Nonrigid Curved Surface • Further changes of surface sound pressure levels are 
expected relative to case (c) where the surface was rigid. 

The main purpose of this brief survey of reflection, or scattering, is to emphasize that the 
sound pressure level perceived on the surface is a function of the impedance perceived by the 
incident field. Quantitative estimates of the effects discussed in Items (c) and (d) above may 
be made using cylindrical scattering analysis. These will be discussed in more detail later and 
extended to include refraction effects. 

23 . 1.12 Refraction 

a. Rigid Flat Surface with an Attached Boundary Layer • In such a case, an incident 
sound field is refracted (there is a change in the direction of propagation relative 
to that of the incident field) in the boundary layer before being reflected by the 
rigid surface as illustrated in Figure 60. Below a certain angle of incidence (the 
so-called critical angle), the incident field may be so refracted that reflection takes 
place in the shear layer above the rigid surface. 

In such an instance, the rigid surface is shielded from the incident field, and diffrac- 
tion analysis must be used to estimate the sound field reaching the rigid surface. 

Such an analysis has been carried out in 2D by Hanson (Reference 2.3.1), McAninch 
(Reference 2.3.4) and Mungur (Reference 2.3.5). Hanson studied the refraction effect for a 
flight configuration using a single and multiple shear layer over a rigid surface. McAninch 
analyzed the refraction in a wind tunnel configuration. Mungur studied both the flight and the 
wind tunnel configurations by numerically integrating a reduced form of the acoustic wave 
equation in the presence of a sheared flow and then matching the incident and reflected fields 
just outside the shear layer. This is equivalent to determining reflection from a surface at the 
edge of the shear layer whose admittance is equal to the numerically transferred admittance 
across the boundray layer of the rigid surface below the boundary layer. 

Figures 61 and 62 summarize results of the 2D analysis, for both the flight and the wind 
tunnel configurations. That significant boundary layer shielding can occur for upstream 
incidence is the main conclusion reached by all three investigators (Hanson, McAninch, and 
Mungur). 
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Figure 60. Convection and Refraction of Incident and Reflected Plane Sound Wave in the Boundary Layer of a 
Rigid Flat Plate in Wind Tunnel Configuration. 



Tunnel Configuration 
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Figure 61. Variation of Wall SPL, 20 log (P w /Pi), Relative to Incident SPL with Angle of Incidence for Mach 
Number of 0.8 and Frequency Parameter ko6 = 1.62 (Jetstar Model Test) and ko6 = 0.469 
(Full Scale). 









Modeling that comes closer to the case of the scattering of propeller noise by a fuselage 
is a 3D scattering problem in the presence of flow and fuselage boundary layer. Hanson 
(Reference 23.6) has studied this case in some detail and shown some significant differences 
relative to the 2D model discussed earlier. This 3D scattering problem is the main subject of 
this work and is discussed in more detail in the next section. 

2.3.2 Scattering of Propeller Noise by a Cylindrical Fuselage 

Consider a propeller of radius Rp and center S (source) mounted on one side of a fuselage 
of radius Rf at a distance Rcl from the fuselage center as shown in Figures 63 and 64. 



Figure 63. Fuselage Scattering of Propeller Noise. 

At an arbitrary field point F, the sound field consists of the sum of a direct field Pd and a 
scattered field Ps. Both the direct and the scattered fields have characteristic directivities 
which are functions of the frequency and other parameters which are defined later. 

Estimates of the amplitude, directivity, and phase of the scattered field relative to an 
incident field. Pine, are made using acoustic diffraction analysis. Geometric ray theory is not 
adequate, due to the wavelength being of the order of the scattering dimension. 

Of equal importance is the sound pressure level distribution on the surface of the fuselage, 
not only in the plane where the propeller is located, but also away from it and all around the 
fuselage, including shadow regions. 

As indicated in Section 23.1, the combined sound pressure level depends not only on the 
amplitude of the incident and of the scattered fields, but also on their relative phase. This rela- 
tive phase is a function of the fuselage surface impedance and the acoustic properties of the 
medium in which the sound is incident and scattered. For a rigid fuselage in a homogeneous, 
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Thickness 


Figure 64. Illustration of Scattering of Propeller Emitted Sound Field by a Cylindrical Fuselage with 
Attached Boundary Layer in Forward Flight. 
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the scattered field on the source side of the fuselage is nearly in phase with the incident field, 
so that the combined field is increased in amplitude. Similarly, as may be expected, in the 
shadow region on the fuselage surface away from the source, the scattered field is nearly in 
antiphase with the incident field. Consequently, the combined field is lower than the incident 
field, thus, producing the so-called shadow region. 

For a rigid fuselage with an attached boundary layer, both the amplitude and the phase of 
the scattered field relative to the incident field will be modified, resulting in a modified 
combined field, relative to the case with no attached boundary layer. In the analysis that 
follows, the boundary layer is treated as part of the fuselage, with a modified acoustic surface 
admittance; in other words, the fuselage surface would not appear rigid to the incident field. 
The estimation of this modified acoustic surface admittance would be part of the scattering 
analysis. 

In this study, reflection, refraction, convection, and diffraction are all part of the scatter- 
ing analysis used. The main objective is to estimate delta -dB (the increase or decrease of the 
pressure field on the surface of the fuselage relative to the free field, that is, the incident field 
that would have been there in the absence of the fuselage). As will be shown, several 
parameters influence the delta -dB. 

2.3.2.1 The Governing Wave Equation 

The sound field emitted from the propeller and incident on the fuselage, and the resulting 
scattered field, must satisfy an appropriate governing wave equation. 

For a mean flow field that is not uniform, such as that in the boundary layer of a cylindrical 
fuselage in flight (Figure 65); the appropriate governing wave equation may be written in the 
form: 


v 2 p + 2p 


3U 3v 
9R 3z 


1 D 2 P 
c 2 Dt 7 


0 


( 2 . 3 . 1 ) 


where 


V 2 = three-dimensional Laplacian in cylindrical coordinates, 
= convected derivative = + U 


U = axial component of the mean flow, which is a function of the radial coordinate R. 
3U 

jjjj = Mean radial shear and is also a function of the radial coordinate. 


v = radial component of the unsteady flow associated with the sound field. 
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A solution for Equation 23.1 may be expressed as: 

= £ A (tu) (R,k ) '-*2 

q q q z 

where Fq (R, k z ) must satisfy a reduced form of Equation 23.1; namely, 


( 2 . 3 . 2 ) 


d 2 F 

w 


i k * S! S ♦{ k .* - V • $\ F * 0 < 2 - 3 - 3 > 


where: 


ko — <ii/Co 
K = kz/ko 
M = U/co 


and 


q = an integer 0, ±1, ±2, . . . 

Each term of Equation 23.2 represents a cylindrical wave spinning around the fuselage, 
propagating downstream or upstream (depending on the sign of the wave-vector k z ), and 
radiating radially away from the surface of the fuselage. 
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In flight, outside the boundary layer, the Mach number, M, becomes zero, and the shear 
term, dM/dR, vanishes as illustrated in Figure 65. In this case, outside the boundary layer, 
Equation 233 simplifies to a more familiar differential wave equation which may be satisfied 
by a Hankel function, so that Fq(Rdtz) in Equation 232 may be replaced by: 

F q (R,k z ) = « q (k R R) (2.3.4) 

in which kR and k z are interrelated by the dispersion equation obtainable from Equation 2.3.3 
in the form: 

k 2 + k 2 = k 2 = k 2 (l-M K) 2 ( 2 . 3 . 5 ) 

R z - o ' o 


By definition, the wave-vector, k, is associated with the direction of propagation of the 
wave-field. The axial and radial components of the wave-vector, k, are given by: 

k = k cos 0 ( 2 . 3 . 6 ) 

z ~ 

and 

k^ = k sin 0. ( 2 . 3 . 7 ) 

where e is the direction of radiation along which the wave-vector k is aligned. 

In flight, as shown in Figure 65, the acoustic medium outside the fuselage boundary layer 
is stationary. From Equation 23.5, with Mo set to zero, 

Ikl = Ik I * Cu»/c ) 

~ ~o o 

However, because the source is moving with Mach number Mo, the frequency perceived 
by a fixed observer outside the fuselage boundary layer is Doppler shifted; thus, 

- * (1+M cos 0) c (2.3.8) 

o o 


hence, 


k = 
z 


k cos 0 = 


(1+M cos 0) 
o 


cos 0 


and 


k R s ^ 


k 

sin 0 = j-. M 0 

(1+M cos 0) 
o 


sin 0 
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It should be noted that the angle e is the angle of emission as shown in Figure 65. It is related 
to the geometric visual angle 0v in the form, 

0 = 0 V - sin -1 (M o sin 0 y ) (2.3.9) 


2.3.2.2 Scattering Analysis 

Having defined the governing wave-equation and the forms that the sound field may take 
outside the boundary layer, we may now refer to Figures 63 and 64 and write down some formal 
solutions. 

The sound field at an arbitrary field point F at location (R,<j>,z,u>), relative to a cylindrical 
coordinate coincident with the fuselage axis, may be written as: 

P(R,$,z,ut) = P d (R,<M,u)) + P s (R,<M,iu) (2.3.10) 


where Pd refers to the direct field emitted from the propeller source, and P s refers to the field 
scattered by the fuselage surface as a result of an incident field Pine. If the field point is on the 
fuselage surface, then the direct field and the incident field are coincident. The next objec- 
tive is to evaluate the scattered field P s in terms of the incident field Pine. The analysis must 
be carried out in the frequency domain, because diffraction and scattering (including refrac- 
tion) are frequency dependent phenomena. 

2332.1 The Scattered Field 

Outside the boundary layer, as discussed in Section 23.2.1, the scattered field may be 
expressed in terms of cylindrical harmonics of yet undetermined amplitudes Aq(o) in the form: 

P s (R,<M,o>) = X A q (©) Hq (k R R) (2.3.11) 


where 

Hq(kRR) = cylindrical Hankel function of integral order and of the second 
kind (outgoing waves) 

= | k I sin 0 
k 2 = |k| cos 0 

k = propagation wave-vector = k o /(l+M Q Cos0) 
k = (u»/c ) 

O 0 
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The amplitudes of Aq of the scattered field as expressed in Equation 23.11 must now be 
evaluated in terms of the incident field Pine as shown in Figure 63. This incident field is that 
emitted from the propeller noise source S at angle e s and is a function of the propeller geometry 
and operating conditions. 

2 3.2.23 The Incident Field 

A spectral component of the sound emitted from a B-bladed propeller may be expressed 
in cylindrical harmonics of known or specifiable amplitudes in the form: 


p e fr'.Be »n) = D„ K,e e 0 (k R R0 e' jW<D »‘ 


(2.3.12) 


where: 

n = harmonic integer of the propeller blade passing frequency 
m = nB 

am ~ ConBMt/Rp 

Co — adiabatic speed of sound 

Rp = propeller radius 

Mt = propeller tip circumferential Mach number 

kR = |k| sin e* 8 

k z = |k| cose'e 

|k| = |ko|/(l + MoCOS 0' 8 ) 

Dn(o>n»0'e) = spectral directivity of the emitted field 
Hm ( ) = cylindrical Hankei function of order m. 

In the above expression for the emitted sound field (r', e' g , 4 >') are spherical polar coordi- 
nates relative to an axis through the propeller center, and (R\ <j>', z') are cylindrical coordinates 
also relative to the propeller center, with z' = r' cos 0' 6 . 

To express the scattered field of Equation 2.3.11 in terms of this incident field and 
boundary conditions on the surface of the fuselage (expressed in terms of a fuselage surface 
admittance), it is more convenient to express the incident field in coordinates coincident with 
the fuselage axis. It may be pointed out that the incident field is the same as the emitted field. 
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Each spectral component of the incident field, as expressed by each term of Equation 
23.12, may be expanded in terms of the (R, z, <j>) coordinates coincident with the fuselage and 
expressed in the following form (Reference 23.7): 


P. 

me 


c v ...) = 


( VV 


e" jk z z H m (l^ R') 


( 2 . 3 . 13 ) 


where: 


= D 


<vV> 


e'j k z 2 


H 


m+£ 


R cl> 






£=-<» 


Jl = an integer 

Rci = the distance between fuselage center line and propeller center 

The relationship between (R, <j>, z) and (R', <t>', z') is shown in Figure 66. 

Equation 23.13 expresses the propeller-emitted sound field in the form of multiple 
cylindrical waves converging onto the cylindrical fuselage. The axial wave number k z of the 
incident wave is the same in either representation, with kz = k cos 6'. 

23333 Boundary Condition at the Fuselage Surface in the Absence of a Boundary Layer 

The scattered field P s in Equations 23.10 or 23.11 associated with each spectral com- 
ponent of the incident field as expressed by Equation 23.13 may now be evaluated by 
subjecting the combined field (incident plus scattered) to pertinent boundary conditions at or 
near the fuselage surface. This may be specified as follows: 

In the absence of a boundary layer and no mean flow: 


at R = Rf, the fuselage radius, 


l c o V R 
P 


where: 


= -AptRy) 


( 2 . 3 . 14 ) 


P = Pine + Ps of Equations 23.13 and 23.11, 

p and c 0 = mean density of the medium and adiabatic speed of sound, respectively, 
v R = radial component of velocity associated with combined sound field P, 
and 

A f(Rf) = fuselage surface acoustic admittance. 
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Figure 66. Cylindrical Coordinate Relative to Fuselage Axis (Forward Looking Aft). 



The radial component of the velocity field, vr, may be expressed in terms of the radial 
derivative of the combined pressure field, using the linearized momentum equation: 


- l dP 

p c o v R ‘ ' jk ffl 
J o 

Substituting Equations 
following: 


(2.3.15) 

23.15, 23.11, and 23.13 into Equation 2.3.14, one may obtain the 


n ^ 


(2.3.16) 


*Z*T \ »„< ) ‘‘ lq * ♦ w ? *F W > J l< ) 

q £ 

where: 

h,< )= ^ {H ’ ( » 


J 7 )= dR W t ( )} 


and the argument in H'q, J^, Hq and J £ is IcrRf, and that in Hm + 1 is IcrRcI. 

From the orthogonality of e‘ jq< ^ and e"^, Equation 2.3.16 may be simplified and 
rearranged in the form: 


Aq H'q ( ) = -5q£D n (o^) ) J j ( ) Yq£ 

where 

{l - )/j- £ ( )} 

V = {1 - AjH q ( )/H' q ( )} 
and 

6 » = Kronecker delta 
q* 

= 1 when q = £ 

= 0 when q * £ . 


(2.3.17) 


(2.3.18) 
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Equation 2.3.17 allows the scattered amplitude, Aq, to be expressed in terms of the 
amplitude D n (on) of the spectral component of the incident field. It may be pointed out that 
if the fuselage surface is rigid, the fuselage surface admittance becomes zero, and Y £ in 
Equation 2.3. 17 becomes unity. q 

In this case, 

Aq = -Dn(can) Hm+q( ) J' q( )/H' q( ) (2.3.19) 

If the fuselage surface is not rigid, then both the amplitude and phase of the scattered 
amplitude are modified relative to that of a rigid fuselage. 

In evaluating the Bessel and Hankel functions, J* ( ) and H q ( ) and their derivatives, it 
may be recalled that the argument common to all is IcrRf. In an ambient medium, this may 
be specified as koRFsinei, where 9i is the angle the incident field makes with respect to the 
propeller axis. 

Equation 2.3.17 will be used to reconstruct Equations 23.10 and 23.11 in terms of the 
incident field amplitude. This part of the analysis is deferred until after a discussion of the 
boundary layer effect. 

2333.4 Effect of the Presence of a Fuselage Boundary Layer on the Boundary 
Conditions 

The boundary condition at the fuselage surface as expressed by Equation 2.3.14 is still 
valid. However, the pressure fields as expressed by Equations 23.11 and 2.3.12 are valid only 
outside the boundary layer and, therefore, cannot be applied at the fuselage surface. 

To overcome this problem, the boundary condition is applied at the outer edge of the 
boundary layer, using an acoustic transferred admittance instead of the actual fuselage surface 
admittance. Such a transferred acoustic admittance is obtained by numerically integrating the 
applicable governing wave-equation; namely Equation 233, from the fuselage surface to the 
edge of the boundary layer. A fourth order Runge-Kutta integration scheme is used. Such a 
scheme requires subdividing the integration region (for our application, the boundary layer 
thickness) into several points, and evaluating the coefficients of the governing equation at and 
between each of these points. Details of such an application are given in Reference 23.8. In 
summary, the pressure field and its radial derivative at the outer edge of the boundary layer 
(Rf + s) may be expressed as: 


I l«M 

where P and P' are the pressure field and its radial derivative at the fuselage surface, respec- 
tively. [T] is a 2 x 2 transfer matrix obtained by the numerical integration of Equation 2.3 .3 
expressed in the form: 


P 

P' 


(2.3.20) 
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(2.3.20b) 


F 

/ 

0 1 


F 

F’ 


-a -b 
• m 


F' 


where ' refers to the derivative with respect to R, and where a and b are the coefficients of F 
and dF/dR of Equation 23.3. 

Evaluation of the transfer functions [T] requires evaluation of the coefficients of Equation 
2 33 from point-to-point across the boundary layer. Consequently, the velocity profile; the 
Mach number, Mo; the spinning integer number, q (0, ±1, ±2, . . .); and (kz/k 0 ), or the angle 
of incidence, must be specified. Thus, Equation 2.3.2 may be evaluated at each specified loca- 
tion on the fuselage relative to the propeller plane. 


At the edge of the boundary layer, the following boundary condition is imposed: 



( 2 . 3 . 21 ) 


as 


e -* 0 

In evaluating (P'/P) at Rf + 8-e, Equation 2.3.20 is used; and in evaluating (P'/P) at 
Rf + 8 + e. Equations 23.11 and 23.13 are used. Thus, from Equation 23.20, 

,,. x *.. ♦ 

‘ P > R f +6-c ‘ T » 7 T 12 (P'/P) R| (2 - 3 - 22 

and (P' /P)Rf is replaced by -jko«4F(RF) from Equations 23.15 and 23.15. 

Equation 2332 is a measure of the fuselage surface transferred admittance and may be 
denoted as -At. 

Using Equation 2.3.2, expressions similar to Equations 2.3.16 and 2.3.17 are obtained with 
At replacing ^1 f. 

The development in this section allows the use of classical scattering analysis to include 
boundary layer refraction effects. 

2323.3 Evaluation of the Sound Pressure Level on the Fuselage Surface 

The sound pressure field on the fluid surface just outside the fuselage boundary layer is 
the sum of the incident and of the scattered fields. Using Equation 23.13 for the incident field 


too 



and Equation 2.3.11 (rearranged somewhat) for the scattered field, the sound pressure level 
at this surface (Rf + 5 , $) may be written in the form, 


- VvV *' ik * z E Vl ( Vci> 3 t M «' J1 * 


( 2 . 3 . 23 ) 


+ L A q H' (a){H q (a)/H' (a)} e ’ jq *“ jk z z 


where 


a = k R (Rf + s), 


and the summations on 1 and q are theoretically from -« to + », but will be truncated due to 
the convergence of the series with increasing order. 


Equation 2.3.17 allows the amplitude and phase of Aq H’q (a) to be expressed in terms of 
the incident amplitude Dn(u>n,e e ) and of the fuselage surface transferred admittance. If this is 
done, then Equation 2.3.23 becomes: 

* VVV e " jk2Z 

( 2 . 3 . 24 ) 


* 



W» 





H»(a) . -I 

V<*> 577) h |. 


where 

b = kRRd 
a = kR(Rf+ 6 ) 

= y ^ of Equation 2.3.18 with q = 1 . 

and k R = jk sine has been discussed previously and is defined in Equation 2.3.7 or 2.3.7a. 

Equation 2.3.24 represents the principal result of the acoustic analysis accounting for the 
modification to an incident field (of the wave emitted from a propeller) by the scattering due 
to a cylindrical fuselage in the presence of convection, refraction, and diffraction. 

Since the result is more conveniently expressed in terms of an increase or decrease relative 
to the incident field, the directivity term D n (wn,0) does not need to be specified; it is canceled 
by the normalization process. 
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To obtain the acoustic pressure field on the surface of the fuselage, Equation 2.3.24 is used 
in conjunction with the pressure component of Equations 2.3.20 or 23.22 relating pressure 
fields at (Rf + 5) and at Rf. Thus, 

P(R f +6,X) ={Tn + T 12 (P'/P) Rf | P(R f ,A) (2.3.25) 

Let P(Rf)/P(Rf + 5) be denoted by PWOBL( l ) for each of the spinning modes i, where 
PWOBLU) s 1/jTu + Ti 2 (P7P) Rf } (2.3.26) 

and is evaluated numerically when integrating Equation 2.3.20b. 

The pressure field on the fuselage surface may now be written as: 

P(R-) = E P(R f ,*) P(R f +6,£) PWOBL(i) (2.3.27) 

1 l 1 l 


where P(Rf+S , l) is the expression in Equation 23.24 inside the sign corresponding to 
s umma tion over %. 

2.3.3 Computed Results 

A computer program has been developed to evaluate the sound pressure field on the 
fuselage surface and on the surface at the outer edge of the fuselage boundary layer, as 
expressed by Equations 2336 and 2337. 

Numerical results have been computed to show the various effects of scattering by a 
fuselage and its attached boundary layer. As mentioned in Section 23.1.1, such scattering 
includes effects of reflection, refraction, and diffraction. These effects are functions of various 
geometric, aerodynamic, and operational parameters. Figures 66 and 67 will facilitate inter- 
pretation of some of the parameters used in presenting the results. 

Unless indicated otherwise, most of the results have been evaluated for an 8-bladed 
propeller rotating at a rotor tip Mach number of 0.8, with a fuselage/propeller radius ratio of 
unity, and zero propeller tip clearance. The fuselage surface is assumed to be rigid. The mean 
velocity profile in die boundary layer corresponds to a Blasius profile on a flat plate. 

2.3.3.1 Reflection and Diffraction in the Absence of Mean Flow 

As a prelude to the presentation of results. Figure 68 has been chosen to show some salient 
features of reflection and diffraction. Two polar plots are presented demonstrating azimuthal 
distributions on and around the fuselage surface of the incident field and of the combined 
incident + scattered field, for sound emitted and incident at 90° to the propeller axis. The 
azimuthal angle is measured at the fuselage center relative to the line joining the propeller 
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Figure 67. Relevant Parameters in the Scattering of Propeller Emitted Sound by a Fuselage. 


center and the fuselage center as shown in the inset and in Figure 66. The incident field is 
normalized to that at <j> = 0, and is actually evaluated as 20 logio { Pine (Rf,<|>)/Pmc (Rf, 0) }. The 
deviation from 0 dB represents the decay of the emitted field from the propeller with distance 
and may differ from the inverse square law due to near field decay. 



Figure 68. Relative Sound Pressure Level Distribution in dB on the Surface of a Rigid Fuselage 
Arising from the Scattering of an Incident Field Emitted from an 8-BIaded Propeller 
Operating at a Tip Mach No. of 0.7. 


The other plot, denoted as incident + scattered fields, is the corresponding SPL distribu- 
tion on the fuselage surface and is relative to the local incident field at the same location. The 
0 dB circle is a reference level. Levels above 0 dB correspond to amplification, or in-phase 
reflection; levels below 0 dB correspond to attenuation or shadowing as a result of out-of- 
phase reflection (scattering). 
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Several features may be observed and commented on: 

• There are two distinct regions on the fuselage surface where the SPL is above and 
below the 0 dB circle. In the region 93° > «|» > -123°, the scattered field construc- 
tively interferes with the incident field to produce levels above 0 dB. In the region 
93° < <j> < -123°, a shadow zone is present where the scattered field interferes 
destructively with the incident field, producing relative levels below 0 dB. 

For this case of fuselage radius equal to propeller radius and zero propeller tip 
clearance, the geometric shadow should be 60° < <|> < -60°. Within the shadow 
region, the SPL varies azimuthally producing “hard” and “soft” shadows which 
correspond to the well-known Fresnel diffraction. 

The depth of the shadow zone and the number of peaks and valleys within it are 
functions of the propeller frequency, blade number, and fuselage radius (not shown 
in the figure). 

• A second important observation is that the fuselage SPL distribution is not 
symmetrical about the line joining the propeller and fuselage centers. The peak 
relative SPL for the case considered is around <t> = -40. This asymmetry is due to 
the fact that the incident field is that emitted from an 8-bladed propeller rotating 
in a counterclockwise direction. This asymmetric feature has also been observed 
by Hanson (Reference 2.3.6). 

• The peak SPL on the fuselage surface at 4> = 0 is less than 6 dB, relative to the inci- 
dent field, supporting the experimental results of Hubbard et al. (Reference 2.3.3) 
that reflection from curved surfaces does not produce pressure doubling. In fact, 
near pressure doubling (5.8 dB instead of 6 dB) occurs at 4 > = 320° = (-40°). 

2.Z.ZJ2 Influence of Fuselage Curvature 

Figure 69 shows the azimuthal distribution of the combined incident and scattered fields 
(normalized by the local incident field) for two fuselage-to-propeller radius ratios, namely 
(Rf/Rp) = 1 and 2. 

The other relevant parameters that are common to both distributions are: 

• Blade number = 8 

• Blade tip Mach number = 0.8 

• No mean flow. 

The main points of observation are: 

• The relative SPL at <j> = 0 increases from 5.0 to 5.4 dB as (RflRp) increases from 
1 to 2 
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Figure 69. Influence of Fuselagc/Propeller Radius Ratio on Combined (Incident and Scattered) 
Relative SPL Distribution on and Around Fuselage Surface for an 8-Bladed Propeller 
with Blade Tip Mach No. of 0J. 





• The peak relative SPL is 5.8 dB at <|> = 335° 

• The shadow region is larger and deeper for the higher radius ratio. 

2.3.3.3 Influence of Angle of Incidence on SPL Azimuthal Distribution 

Figure 70 shows the no mean flow relative SPL distribution (incident + scattered fields 
nor maliz ed to incident field) on and around the fuselage surface for different angles of 
incidence. 

Points of observation are: 

• The azimuthal distribution is not symmetrical 

• The shadow region reduces for lower angles of incidence 

• For 9j £ 50°, regions of amplification exceed 6 dB. 

An amplification region (or regional focusing) prior to a shadow region is not uncommon 
in diffraction of wave fields such as that from a knife edge. Such a feature is apparent 
for 0j = 50°; however, the degree of amplification is unexpected. As will be shown later, the 
presence of flow and boundary layer modify these distributions. 

2.3.3.4 Influence of Mean Row and Boundary Layer 

An example of the influence of mean flow and boundary layer on the azimuthal distri- 
bution of the combined incident + scattered fields (normalized by the incident field) is shown 
in Figure 71. It corresponds to a flight configuration of Mach No. 0.8 with angles of incidence 
of 50° and 130°. The boundary layer thickness is 0.01 of the fuselage radius. The blade number 
and the tip rotational Mach number are 8 and 0.8, respectively. The corresponding distribution 
for zero cruise Mach number is also shown. 

All three curves correspond to azimuthal distributions around the fuselage at the outer 
edge of the boundary layer. The distribution at the fuselage surface is the same as that at the 
outer edge for the zero flight Mach number and for M c = 0.8, with 0i = 130°. For 0i = 50° 
and M c = 0.8, the fuselage surface pressure at R = Rf is 02 dB lower than that at R = Rf + 6 
(the outer edge of the fuselage boundary layer). This difference is larger for thicker boundary 
layers as is shown later. 

The main points of observation are: 

• Flight Mach number alters the symmetry of the combined pressure distribution 
relative to angle of incidence 

• The flight Mach number reduces the relative SPL at <|> = 0 from 5.2 dB to 4.1 dB 
for 0i = 50°, and to 4.8 dB for 0i * 130° 
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Figure 70. Effect of Angle-oMncidei.ee on the Relative Combined Incident and Scattered SPL 
Distribution on and Around the Fuselage Surface. 
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• Slight shadowing effect, -0.2 dB (not shown) for SPL below fuselage boundary 
layer. 

2.3.3.5 Influence of Fuselage Boundary Layer Thickness and Cruise 
Mach Number 

The influence of refraction due to the presence of a fuselage boundary layer is shown in 
Figure 72. The combined sound pressure level (incident + scattered fields normalized by the 
incident field) at the azimuthal angle <|> » 0 is shown as a function of visual angle of incidence. 
Two plots are shown, corresponding to levels above and below the fuselage boundary layer; 
namely at R = Rf + 5 and R = Rf. The computations have been done for a flight Mach num- 
ber of 0.8 and a boundary layer thickness ratio of s/Rf = 0.10. The actual profile across the 
boundary layer is that corresponding to a Blasius profile on a flat plate. 

The main points of observation are: 

• SPL relative to free field above the fuselage boundary layer and at <j> = 0 is higher 
than 0 dB (amplification) for all angles of incidence, ranging from 1.0 to 4.8 dB 

• SPL relative to free field on the fuselage surface below the boundary layer is lower 
for upstream incidence than for downstream incidence and corresponds to bound- 
ary layer shielding 

• Below the fuselage boundary layer, the range of angles of incidence where 
amplification (level above free field) exists is limited between 92° < 9i < 140°. 

Figure 73 shows the relative SPL distributions on the fuselage surface for three different 
boundary layer thickness ratios; namely S/Rf = 0.01, 0.05, and 0.1. The trend of increased 
refraction as (s/Rf) increases is as expected. 

Figure 74 illustrates the influence of cruise Mach number on the relative SPL distribution 
at $ = 0 ad on the fuselage surface for s/Rf » 0.1. Again, the trend of increased fuselage 
boundary layer shielding is as expected. 

The case of cruise Mach number equal to 0.8 compares favorably with Hanson’s analytical 
results (Reference 23.6), allowing for possible differences in boundary layer profile and 
boundary layer thickness. It may be noted, however, that the other two cases (corresponding 
to cruise Mach numbers of 0.6 and 0.7) do not show the trends indicated by Hanson; namely, 
no refraction beyond 2-feet forward of the propeller plane. The present results show consis- 
tent trends of shielding similar to the case of Mach No. 0.8, but with less refraction. 

2.3.3.6 Acoustic Shielding at Higher Harmonics of Blade Passing Frequency 

All of the results presented thus far have been for the case of one times blade passing 
frequency. Figure 75 shows similar relative SPL distributions for three harmonics of the blade 
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Figure 72. SPL Distribution (Incident and Scattered) Relative to Incident Field Below and Above 
Boundary Layer. 


RELATIVE DB: 20Log{(P. +PJ/P. 



Figure 73. Influence of Boundary Layer Thickness on Relative SPL Distribution on 
Fuselage Surface at $ 3 0. 
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Figure 75. r>nhln«d Incident and Scattered Relative Sound Pressure Level Distribution on 
the Fuselage Surface at <J> = 0 for First Three Harmonics of BPF. 
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passing frequency. These distributions are the combined SPL (incident + scattered field 
normalized by the incident field) on the fuselage surface at R = Rf and <t> = 0. 

As expected, the combined fields at the higher harmonics undergo more refractive shield- 
ing. An interesting observation is that for angles of incidence 9i > 90°, the combined sound 
field at the blade passing frequency undergoes very little shielding; this is due to the thin 
boundary layer chosen for the illustration, namely s/Rf = 0.01. The corresponding fields at 
the higher harmonics undergo more shielding. 

2.3.3.7 Comparison of 2D Plane Wave and 3D Cylindrical Wave 
Scattering by Boundary Layer 

In Subsection 23.1.12, the refraction of a 2D plane wave by a rigid flat surface with an 
attached boundary layer was discussed. Figure 62 shows the boundary layer shielding effect 
for two frequencies. 

To make a comparison with this 2D analysis, the cylindrical scattering analysis has been 
used to compute the SPL distribution for the following case: 

Number of blades (B) = 8 

Rotational tip Mach number (Mt) = 0.8 


Cruise speed M c = M z = 0.8 
Fuselage/propeller radius ratio (Rfi'Rp) = 1 
Boundary layer/fuselage radius ratio s/Rf = 0.073 
The wave number: 

R f A 

k 6 = k R Jj- 
o o p R R 


* BN t x 1 x 0.073 


= 0.467 

This wave number contrasts with that of 0.469 for the 2D analysis. 

Figure 76 portrays the relative SPL distribution at <j> = 0 and R = Rf. 

The 2D plane results are also presented; demonstrating remarkably similar trends. It 
appears that the 3D analysis yields lower amplification than the 2D case. In the shielding 
region, 9 < 80, the 3D analysis yields somewhat more shielding than the 2D analysis. It must 
be stressed that this comparison is for <j> = 0. As pointed out in Subsection 2.3.3. 1, the 
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amplification is higher at some other value of 4>. However, for <(> * 0, comparison with 2D 
plane wave refraction loses its meaning. 

2.3.4 Concluding Remarks on Fuselage Scattering 

An acoustic scattering analysis has been developed to make estimates of the acoustic 
shielding and amplification effects associated with sound incident upon a cylindrical fuselage 
supporting a mean flow boundary layer. 

The es timat e is evaluated in the form of a correction to be added to or subtracted from 
the incident free field. The incident field is chosen to be representative of acoustic emission 
from a propeller of specifiable blade number and rotating tip Mach number. 

The main conclusions are as follows: 

1. The combined sound pressure level distributions (incident + scattered) on the 
fuselage surface and at the outer edge of the boundary layer are not azimuthally 
symmetrical about the line joining the propeller and the fuselage center lines. 

2. Regions of relative amplification and relative attenuation, corresponding to 
ill umina ted and shadow regions are functions of the fuselage-to-propeller radius 
ratio, the mean flow or cruise Mach number, the boundary layer thickness, and the 
harmonic of blade passing frequency under consideration. 

3. On the fuselage surface, acoustic shielding occurs even in the illuminated region 
(<j> = 0), especially for angles of incidence 0i < 90°. 

4. This refractive shieldin g becomes larger with increasing boundary layer thickness, 
mean flow or cruise Mach number, or frequency of the sound field. 
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3.0 INSTALLED TURBOPROP NOISE ANALYSIS 


3.1 Configuration 

3.1.1 Background 

Under a separate NASA-sponsored effort, the PTA (Propfan Test Assessment) Program, 
acoustic data have been taken (under various flight conditions) from a large-scale propfan, 
mounted on the wing of a Gulfstream aircraft. Although data/theory comparisons were not 
required for the installed turboprop noise analysis portion of the current contract, it was felt 
that exercising the prediction models developed in the course of this work in order to predict 
the noise generated by a real (rather than a fictitious) propfan/airframe combination would 
be more me aningf ul. Consequently, it was decided that predictions should be made under 
conditions corresponding to those of three of the flights that were made during the test 
program, and that the predictions should be made at locations corresponding to those at which 
measurements had been taken on the aircraft. The three flight conditions selected represent 
the design point for the propfan (high altitude, high speed), together with two low altitude, 
low flight Mach number points at differing power and tip speed. 

3.1.2 The PTA Airplane 

The large-scale, 9-foot-diameter, SR-7L propfan was mounted on the left wing of a 
modified Gulfstream G-II business jet transport and was flight tested by the Georgia Division 
of the Lockheed Aeronautical Systems Company, under the NASA-sponsored Propfan Test 
Assessment Program (Reference 3.1.1). During this test, acoustic data were acquired with 
microphones located both on the fuselage of the airplane and on a boom situated further out- 
board than the propfan, and at the same sideline distance as the closest points on the fuselage. 
Figures 77 and 78 depict the PTA aircraft front view and plan view, respectively, while the 
various microphone locations on the fuselage surface and on the boom are shown in Figures 
79 and 80. Tables 2 and 3 identify the microphone locations, in terms of the airplane coor- 
dinate system, as functions of: fuselage station (axial), water line (vertical), and butt line 
(horizontal). In addition. Table 4 shows the location of the center of the propfan disk in the 
same coordinates. It will be noted that Table 4 has three sets of coordinates, corresponding 
to three values of “nacelle tilt angle.” This arises from the manner in which the propfan was 
mounted on the aircraft. One of the variables investigated during the flight test was the flow 
inflow angle in the vertical plane with regard to the propfan disk; thus, varying the tilt of the 
nacelle in the vertical plane permitted a greater range of inflow angle over the flight envelope 
than would have been possible otherwise. 

The computer programs used to calculate the noise of the installed propfan require 
sideline distance and observer (or emission) angles relative to the center of the propfan disk 
as input, rather than airplane coordinates. Table 5 shows the locations of the microphones on 
the left side of the fuselage and on the boom in terms of sideline distances (DIST in Table 5) 
and observer angles (THETA) from the center of the disk at the two values of nacelle tilt angle 
for which the calculations were required. Also shown are tabulations of PHI, the azimuthal 
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Figure 77. PTA Aircraft - Front View. 
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Figure 78. PTA Aircraft - Plan View. 


M24- 




(2) MICROPHONE DIAPHRAGMS ARE 
VERTICAL AND ON PROPELLER 
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Figure 80. Acoustic Boom Installation and Microphone Locations. 



Table 2. Fuselage Surface Locations for Prediction Purposes • 45 Locations. 
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Table 3. Wing Boom Surface Locations for Prediction Purposes. 


FUSELAGE STATION 
t AXIAL) 

247 

301 

328 

355 

409 



40. 5D 

PROP 

PLANS 

-0.25D 

-0.50 

-1.00 


HATER LINE 

283.965 

283.965 

283.965 

283.965 

283.965 

1 

1 

MJTT LINE 
LOCATIONS 
TOR BOOH 


Table 4. Propeller Disk Center Coordinates for the 
Three Nacelle Tilt Angles. 


NACELLE 

TILT 

ANGLE 

BUTT 

LINE 

BL 

FUSELAGE 

STATION 

FS 

WATER 

LINE 

WL 

+2* 

165 

.3L 

300.235 

77.718 

-1* 



301.478 

72.603 

-3* 


1 

302.454 

69.233 


I 


124 
























Table 5. Microphone Locations in the Coordinates Required by the Computer Programs. 
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angle, in the coordinate systems required by the installation effects and fuselage scattering 
computer codes. 

3 J2 Calculations 

3.2.1 The Cases Analyzed 

Tables 6 through 8 provide details of the propfan operating conditions for each case 
analyzed, and Table 9 defines the abbreviations utilized in these tables. 

Case No. 1 (Table 6) corresponds to the design point for the SR-7L propfan; 35000-feet 
altitude, at a flight Mach number of 0.8, and a tip rotational Mach number also of 0.8. Cases 
2 and 3 are low altitude flights; nominally 1000-feet altitude, and a flight Mach number of 0.3. 
Case 2 has a high tip rotational Mach number of 0.73, with a high power absorption of 5886 
shp. On the other hand, Case 3 has a tip rotational Mach number of 0.62, with a corres- 
pondingly lower power, 4728 shp. 

For each case, both uninstall ed (free field, no unsteady loads included in the calculation) 
and installed (free field, with unsteady loads resulting from: angie-of-attack operation, the 
presence of the fuselage, and the effect of the wing lifting line) predictions were made. The 
uninstalle d cal cu l ati ons were made for the microphones on the boom, as well as on the 93.65 
and 91.71 water lines along the fuselage. The installed calculations were made for the 
microphones around the fuselage in the plane of the propfan (Fuselage Station 301) and at 
0.25 propfan diameters upstream and downstream (Fuselage Stations 274 and 328). In addi- 
tion, for Case 1 (the high flight Mach number case), fuselage scattering predictions were made 
for three fuselage boundary layer thicknesses. 

3.2.2 Modifications to the Codes 

The SRPFAN (Single Rotation PropFAn Noise) prediction code delivered under this 
contract requires the spanwise load distribution on the blades to be input as a swirl (or work) 
coefficient distribution, as described in Subsection 2.1.23 of this report These distributions, 
along with the chordwise loading parameters, are input to the code by means of unique Load- 
files, obtained from an Euler flow solution of the particular operating conditions under 
consideration. For the particular cases reported here, NASA Lewis Research Center kindly 
supplied GE Aircraft Engines with Loadfiles generated from a Denton Euler code (Reference 
32 . 1 ). 

However, during the interim between the delivery of SRPFAN to NASA and the start of 
the present calculations, modifications had been made to the NASA version of the code, such 
that it now required spanwise loading to be input as a distribution of lift coefficient. 
Consequently the Loadfiles received contained lift (rather than swirl) coefficients. The GEAE 
code was modified to accept lift coefficient as input, while retaining the ability to scale the 
input distribution if necessary to match the required power input 
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Table 6. Prediction Case No. 1. 
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( 91 ) PSTQ 


9547.14 

FT-LSF 

31 JTTHH 

• *14.9314 

DCGSES 

e 

{ 94 ) FHST 


1543.14 

LSF 

1 2 ) TTS 0 

• 0.114700 

016111 

c 

( 99) PSMP 


344.040 

HF 

33 )OATC 

• *47.4200 

016111 

6 

( 90PPRT 


1193.39 

FT-LSF 

34) PA 

• 499.023 

FSP 


< 97 ) PPBP 


197.493 

LSF 

39) PTOT 

• 741.242 

FSF 


( 99 ) POPL 


37.9324 

HP/FT* 2 

30 COO 

• 947.734 

FFS 


( 99IJP 


3.12131 


37) SOMA 

• 0.301300 



(100ICP 


1.74714 


34IPHO 

• 0.7144301*03 

SLUOS/FT* 3 

(lOl)CTP 


0.424520 


39IPNOC 

• 0.707437 

LSF-SSC/FT* 3 

( 102 ) PPM 


0.741440 

% 

40 )QR 

• 223.944 

FSF 


(10DRPP1 


223.349 

HI 

4DALPH 

• 1.40954 

016111 





> 

4 2 ) ALMX 

• 2.34490 

016111 





4* 

4] ) AUDI 

• 1.37940 

016111 



A. 



44 ) ALSO 

• 0.220440 

016111 





y 

43 ) 33NB 

• *.472310 

oiaiss 


(109)RPM1 


44.5344 

% 

40S3MX 

• 0.204100 

016111 


( 109 ) R1MX 


44.3304 

% 

47 ) SSMH 

• *.933400 

016111 


(110IR1MN 


44.5000 

% 

40 3330 

• 0.209730 

016111 


(111 1 R130 


0 .1204001*01 

% 

49 INTA 

• *1.00000 

016111 


(112ISPR1 


2072.49 

LSF 

30 ) DKOT 

• *.139720 

FFS 


(113 )LSP1 


0.91472 4 S + 04 

HI 

3 1 ) PCAO 

• 12.0030 

FSZ 


(U4ILST1 


0. 3421041+09 

HZ 

32 ) PCHX 

• 12.0049 

FSZ 


1 1 1 n nrrn 



y 

S3IPCMM 

• 11.9941 

FSZ 



V 


y 

34 ) PCSB 

• 0.1430001*02 

FSZ 


♦ 1 1"* \ iunw 



^ DFR 

S3 I COP 

• 4.33933 

FSZ 





y 

SOTCAR 

• 20.4730 

016111 

c 

(119IRPH2 


92.4955 

% 

37 ) CC 

* 1127.32 

FFS 


(120) R2MX 


92.9334 

% 

so isxac 

• 0.0000001+00 



(121IR2HR 


92.4433 

% 

39) PC 

• 0.1903241*02 

SLU6S/FT* 3 

(122) R2SO 


0.1433001*01 

% 

oo)PMce 

• 2.14420 

LIP* S SC/ FT* 3 

(123 I5PH2 


3043.23 

LSF 

0DTQK3 

• 1430.32 

FT-l.Br 


( 124 IRSP1 


0 . 9 L 4 7 2 31 + 04 

HZ 

(2)T0MX 

- 1433.30 

FT-LSF 


( 129IRST1 


0 . 3421041+09 

HZ 

03 ITQMH 

• 1425.23 

FT-LSF 


(129) • 0. 

00000014*00 0 
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Table 7. Prediction Case No. 2. 


PIT 43. RUM 03R2 13-HOV-47 



1 ) OAT* 

• 1009.47 


( 64 

)TQSD 

• 24.1794 

FT- LSP 


2 ) plt 

- 49.0000 


< 65 

) PSMP 

- 5953.11 

RP 


3) HUM 

• 9.14000 


< 66 

) SNPC 

• 0 • OOOOOOE+OO 

% 


4} STM* 

• 04:31:17.401 

R/R/s 

( 47 

) RPRL 

• 100.256 

% 


9 ) CTHS 

• 04:32:00.401 

R/R/s 

( 64 

) RLM 

• 100.746 

% 


4)TAM* 

• 4.40000 

ocajtcc c 

( 69 

) RLMR 

« 100.070 

% 


7) TAMS 

• 14.9700 

P9Z 

( 70 

) RLSD 

- 0.116900 

% 


4>R* 

- 49.0000 

% 

( 71 

) PNCP 


% 


9 IVSJIX 

- 4.00000 

MOT 

( 72 

) RPT 

• 11529.4 

RPR 


10)W3HR 

• 0.400000 

MOT 

( 73 

) RPNI 

■ 96.3429 

% 


inwsAa 

• 4.30000 

MOT 

( 74 

) RUHR 

- 94.4570 

% 


1 2 ) WS D 

« 399.000 

DEGRSE 

( 79 

) RUHR 

• 99.9444 

% 


13) f AMO 

• 0 . 247749*— 02 

slk/gs/pt*3 

( 76 

} RRSD 

• 0.41 4200**01 

% 


l4jca 

• 1104.10 

PPS 

( 77 

) RPNC 

• 13779.9 

RPR . 


15>PRCG 

• 2.73994 

LSP-S*C/PT*3 

< 74 

)CP1 

- 4363.63 

HE 


14 ) AC GW 

• 41329.0 

LIP 

( 79 

)TP1 

- 13435.3 

HZ 


17)RTRU 

■ 0.303400 


( 40 

) PRPH 

• 1696.16 

RPR 


14 ) HAUL 

• 1030.00 

PC ST 

( 41 

) TROT 

• 799.295 

PPS 


i9)*rr 

■ 547.544 

rs*T 

( 42 

JVPWD 

• 334.095 

PPS 


20)HPM 

- 479.424 

PSST 

( 43 

) TTH 

- 466.309 

PPS 


21 ) HfTH 

• 334.449 

PCST 

( 44 

) PTHR 

• 0.726247 



22 ) *930 

• 149.447 

FCST 

( 45 

)nr*x> 

• 0.303400 



2 3 ) XT AS 

• 197.404 

SHOT 

< 46 

)TMTH 

• 0.747141 



24 IKTM 

• 200.444 

MOT 

( 47 

) SAPP 

- 36.4441 

DEGREE 


29 ) KTRR 

• 192.790 

MOT 

< 44 

) BAM 

- 37.1094 

DEGREE 


24 ) XTSO 

- 2.20074 

MOT 

< 49 

)SAJUI 

• 35.4046 

OCGRCX 


27) TOO 

• 334.099 

PP* 

( 90 

) SABO 

• 0.244670 

DEGREE 


2 4 ) KCAS 

• 194.414 

MOT 

( 91 

)XA 

• 3.70003 

DEGREE 


2 9 ) TTH* 

• 11.9240 

ocamss c 

( 92 

) SHPA 

• 5446.03 

HP 


30 )TTM 

• 12.1719 

OEORCS C 

( 93 

) PSTQ 

• 14225.9 

PT-LSP 


3DTTJW 

- 11.4922 

DCQI.SC C 

( 94 

) PRST 

- 6612.03 

LSP 


3 2 ) TTSO 

• 0.147190 

DSOSES c 

( 95 

) PBRP 

• 735.754 

HP 


33JOATC 

• 4.77444 

dcgrex c 

< 94 

) PPST 

- 2274.24 

PT-LSP 


34) PA 

• 2071.49 

PSP 

( 97 

) PPBP 

- 426.503 

LSP 


3 9 ) PTOT 

- 2204.29 

PSP 

( 94 

) PDPL 

• 72.6471 

HP/FT* 2 


34) COO 
37 ) SOMA 

• 1100.52 

• 1.00774 

PPS 

( 99 
(100 

) JP 
)CP 

• 1.31304 

• 1.01303 



3 4 ) PHO 

• 0 . 239443**02 

SLUOS/PT*! 

(101 

) CTP 

• 0.526410 



39 ) PHOC 

- 2.43759 

LSP— SEC/PT * 3 

(102 

) PPR 

• 44.2324 

% 


40 )Q* 

- 133.494 

PSP 

(103 

) BPP1 

• 226.154 

HE 


41JALPR 

• 4.70003 

DEGREE 

(104 

) RPM3 

• 35.4340 

% 


42) ALM 

• 5.47040 

DEORCX 

(105 

)R3M 

• 35.9219 

% 


43 )AUUV 

• 3.27420 

DEGREE 

(104 

) R3MR 

• 35.7461 

% 


44 ) ALIO 

• 0.444430 

OKGRES 

(107 

) R3SD 

• 0.374000C-01 

% 


4 9 ) 3S1VB 

- -2.35071 

DEGREE 

(104 

)RPH1 

• 64.5379 

% 


44 ) SSR3C 

- -1.47990 

DEGREC 

(109 

) RIM 

- 64.5977 

% 


47 ) ssniv 

• -4.74440 

DEGREE 

(110 

) RlRR 

• 64.4443 

% 


44)3330 

• 0.774740 

DEGREX 

(111 

) R1SD 

• 0 . 27 400 OE-0 1 

% 


49 ) RTA 

• -1.00000 

OEORCS 

(112 

) SPR1 

• 253.452 

LSP 


90 ) OH OT 

• 11.2407 

PPS 

(113 

)LSP1 

- 1002.52 

HZ 


51 )PCA» 

• 14.7940 

PSZ 

(114 

) LST1 

• 3959.94 

HE 


52) pen* 

• 14.1144 

PSZ 

(115 

) RPH4 

• 32.4520 

% 


53 )PCH2V 

• 14.7725. 

PSZ 

(116 

)R4M 

• 32.5352 

% 


54 ) PCSO 

• 0 . $10000*— 02 

PSZ 

(117 

) R4MR 

- 32.3516 

% 


55 ) COP 

• 0.407255 

PSZ 

(114 

) R4SO 

• 0.474300C— 01 

% 


54 ) TCAS 

• 14.4770 

OEORCS C 

(119 

) RPH2 

• 60.2445 

% 


57 ) CC 

• 1119.44 * 

PPS 

(120 

)R2M 

- 60.3242 

% 


54 ) 3XOC 

- 1.00144 


(121 

) R2HR 

• 60.1444 

% 


59) PC 

• 0 • 234 190C-02 

SLUGS/FT* 3 

(122 

) R2SD 

• 0.42 4 900C-0 1 

% 


40 ) PRCC 

- 2.64434 

LSF-SEC/PT* 3 

(123 

) SPR2 

• 154.243 

LSP 


6UTQE3 

- 2711.49 

PT-LBP 

(124 

) RSP1 

- 907.499 

HE 


42) TQM 

• 2724.25 

PT-LRP 

(125 

) RST1 

> 3546.20 

HZ 


43 )TQHR 

• 2592.00 

FT- LSP 

(126 

) RZCQ 

• 1.00594 

G 
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Table 8. Prediction Case No. 3. 


pit <1. awn l» io-hov-«7 



1 ) OATS 

• 

909.4*9 


( <4)TQSD 

• 4.29611 

FT-LSF 


2 ) PLT 

m 

41.0000 


( 63IP3HP 

• 4776.30 

HF 


Dawn 

m 

13.1400 


( CSISNPC 

• 44444.0 

% 


4 ) STMS 

m 

04:34:24.400 

n/n/9 

( <7)APMZ. 

• 47.4554 

% 


S ) STMS 

m 

00:37:14.401 

h/n/s 

( atiautx 

• 44.1436 

« 


SITAMS 

m 

20.4000 

scones e 

( (9 ) autn 

• 47.6914 

% 


7 ) PAHS 

m 

14.4300 

PSX 

( 70) ALSO 

- 0.7672000*41 

% 


Dan 

m 

77.0000 

% 

( 71IPHCP 

• 43.5744 

% 


*)WSMX 

m 


ROOTS 

< 72 )BPT 

• 10103.4 

KPN 


io matin 

m 


KNOTS 

( 73 )aPMH 

• 96.6103 

% 


ii ms ao 

m 


KNOTS 

( 74 ) 1HMX 

• 96.6992 

% 


12 mao 

m 


OCOSOS 

( 7S)BKMH 

• 94.4531 

% 


1 3 ) a AMD 

m 

0.2334140*02 

SLUOS/FT* 3 

( 7S)anSD 

• 0.3949000*01 

% 


14) CO 

m 

1127.30 

PFI 

< 77)«PMC 

• 13413.3 

KPN 


ISIPHCQ 

m 

2.43404 

lsf-scc/ft # 3 

( 7S)CP1 

• 3199.40 

HZ 


ISIACQV 

m 

40137.9 

UP 

( 79ITP1 

• 12124.0 

HZ 


i7)«rao 

m 

0.303290 


( I0IPKPM 

• 1446.36 

KPN 


10 ) NAOL 

m 

990.000 

FOOT 

< ai)vaoT 

• 700.431 

FFS 


19 )HPT 

m 

409.441 

FOOT 

( «2)VPWO 

» 342.747 

rps 


ZOINPMX 

m 

1034.23 

FOOT 

( 1 3 ) VTH 

• 779.794 

FFS 


2DMPTH 

m 

432.237 

FOOT 

( 14) PTH A 

• 0.624233 



22)HPSD 

m 

112.774 

FOOT 

< asiMpno 

• 0.305290 



23) KTAS 

m 

202.929 

KNOT 

( aOTMTN 

• 0.694963 



24)aTNX 

m 

204.932 

KNOT 

( *7) SAPP 

• 41.7465 

OOQKOO 


2 S ) STNH 

m 

200.147 

KNOT 

( aaiSAHX 

« 42.0504 

DEGREE 


20 ) KTSO 

m 

1.17919 

KNOTS 

( apiaAtw 

• 41.3944 

DEGREE 


27) TOO 

m 

342.747 

FFS 

( 99 ) BASO 

- 0.149100 

DEGREE 


2SISCAS 

m 

199.047 

KNOT 

( 91 ) ZA 

• 5.94331 

oeascc 


29 )TTnS 

m 

23.2430 

000000 c 

( 92 ) SNPA 

• 4724.73 

HF 


10 ITTHX 

m 

23.3400 

ocooss e 

( 93 ) PSTQ 

• 16709.1 

FT-LSF 


3 1 ) TTMH 

m 

22.7434 

ocooco c 

( 94) PHST 

• 5191.54 

LSF 


32 )TTSD 

m 

0.243210 

ocoooo e 

( 9SIPBHP 

- 391.091 

HF 


33) OATC 

m 

17.4302 

ocooco c 

( 90PPST 

• 2044.64 

FT-LSF 


34) PA 

m 

2033.04 

FSF 

( 97 ) PPSP 

• 644.947 

LSF 


1 3 ) PTOT 

on 

2192.29 

FSF 

( 9SI9DPL 

• 34.3794 

HP/FT* 2 


30 COO 
17) SOMA 

m 

m 

1122.07 

0.941440 

FFS 

( 991JP 

(100ICP 

• 1.53721 

• 1.26730 



30 ) PMO 

m 

0.2244970*02 

S|pUOS/FT*3 

(101 HTTP 

• 0.364040 



39)PHOe 

m 

2.34413 

LSF*S0C/FT # 3 

( 102 ) PPH 

• 0.644169 

% 


40 )QS 

m 

134.072 

FSF 

( 103 ) BPP1 

• 194.142 

HZ 


41IALPH 

m 

3.94331 

OCOOCO 

( 104 ) APM3 

• 34.9100 

% 


4 2 ) AUfX 

an 

3.24900 

oooocc 

( 109 ) A3MX 

• 0.0000000+00 

% 


43 IAUUI 

on 

1.44070 

ocooco 

(104IA3HH 

• o . ooooooc+oo 

% 


44 ) ALSO 

• 

0.314130 

ocoooo 

( 107 ) B3SO 

• 0.0000000+00 

% 


49)ssna 

■ 

0.2440000*02 

ocooco 

(lOS)APHl 

• 63.7793 

% 


40ISSHX 

m 

0.424000 

ocooco 

( 109 ) AIMS 

• 43.4672 

% 


47)SSHn 

m 

*.334600 

ocoooo 

(110) AIMS 

• 63.7266 

% 


40 ISSSO 

m 

0.264240 

ocooco 

(111 IA1SD 

• 0 . 2740000*01 

% 


49 IOTA 

m 

2.00000 

ocooco 

( 112 )SPH1 

• 207.190 

LSF 


so i oner 

m 

7.34104 

FFS 

(113M.SP1 

• 97666.5 

HZ 


SI IPCAS 

m 

14.4030 

PSX 

( 114 ) LST1 

• 343743. 

HZ 


521PCMX 

m 

14.4271 

FSX 

(119IAPM4 

• 31.2420 

% 


S3 ) PCMH 

m 

14.7432 

FSX 

( 11( ) A4MX 

• 31.3004 

% 


34 ) PCSO 

m 

0.4320000*02 

FSX 

( 117)X4HM 

• 31.1640 

% 


39) COP 

m 

0.331740 

FSX 

(11SIA4SO 

• 0.21 29000*0 1 

% 


30 ITeAS 

m 

22.1340 

ocooco c 

( 119 ) APM2 

• 54.9542 

% 


S7)ce 

on 

1130.36 

FFS 

(129 ) A2MX 

• 59.2300 

% 


sotszoe 

• 

0. 0000000+00 


(121 )A2MH 

« 49.1094 

« 


99) PC 

m 

0.2337490*02 

SLUGS/FT* 3 

( 122 ) A2SD 

• 0.139650 

% 


00 ) PMCC 

an 

2.44221 

LSF— SEC/FT* 3 

( 123 ) SPH2 

• 113.625 

LSF 


C1ITQS3 

m 

2443.01 

FT-LSF 

(124 )ASP1 

• 47404.7 

HZ 


02 ) TQNX 

m 

2492.00 

FT*lpSF 

( 129 I ASTI 

• 343249. 

HZ 


0 3 ) TQNN 

at 

2466.23 

FT*LBF 

( 12* (tlZCO 

• 1.01737 

G 


129 



I 

Table 9. Aircraft and Propeller Operating Conditions. 


Definition of Abbreviations 


F 


s 

5 

u 

o 


A3 

U 

o 

u 


-Q 

CO 

U 


1- OATE-TEST OATS 

2- FLT •PLIGHT NUMBER 

3- RUN •HUH CVUMSCH 

4- STHC-FROM T04 PILE 
j»tTHt»PHOW T04 PILE — _ 

4-TAMBaGROUNO AMBIENT TEMPERATURE 
7-PAMBaGROUNO AMBIENT PHCSSUHC 
B-RM •GHOUHO RELATIVE HUMIDITY 
9-WSMX-QROUNO WIND SPEED NAS 
O-WSMN-GROUNO HI HO iUIO MIN 
X-WSAG-GROUND WIND SKID AVERAGE 

2- WSO •GHOUHO WINO SPEED OIHKCTZOH 

3- PAM DaG ROUND AZH OCHS ITT 

4- CO •GHOUHO SPEED Of SOUHO 

5- PMCOaQROUND CHAH. IMPEOAHCE 
^-TA^dW«A!*C*A/f GHOSH WEIGHT ' 
7-MTRU-FREE3TREAH MACH 

B— MAG L— ALTITUDE AMOVE GHOUHO 
• -HPT •ALTITUDE ABV SEA LEVEL ME AH 
0—HPMX« ALTITUDE ABOVE SEA LEVEL-MAX 
X-MPTMa ALTITUDE ABOVE SEA LSVEL-HH 

2 - M PSD* ALTITUDE ABOVE SEA LEVEL-SO 

3- XTAS-TRUS AIH SPEED-MEAN 

4- KTMX-TRUE AZH SPEED-MAX 
3— RTHRaTRUE AZH SPEEO-HZH 
C-KTSOaTRUS AZH SPESO-SO 
7-VOO •THUE AIH SPEED 
B-KCAS-CALIBRATEO AIH SPEED 

9 -TTNB- FREE STREAM TOTAL TEMP— HEAD 
0— TTMX—PREESTRSAfl TOTAL TEMP-MAX 
X-TTMN«FRCE3TRCAM TOTAL TEMP-MXH 

2- TTSOaFREESTREAM TOTAL TEMP-BO 

3- OATCaFRCESTREAJt STATIC TEMP 

4- PA •PHEESTHEAM STATIC PHESSUHE 

5- PTOT«PHEESTHEAM TOTAL PHESSUHE 
4-COO •PHEESTHEAM SPEED OP SOUHO 
7-SGHA«0CHSITY HATXO 

B-PMO aPHEESTHEAM AZH OEHSZTT 
f-PMOC* PHEESTHEAM CHAH. tMPEOAHCS 

0- GE aOYHAMZC PHESSUHE 

1- ALPNaAMGLE OP ATTACH— HEAD 
Z-ALMXaAHGLS OP ATTACK-MAX 

3— ALMHaAHGLE OP ATTACX-MXH 

4 — ALSO* ANGLE OP ATTACX-SO 

5— SSNSaSIOESLZP ANGLE-MEAN 
B-SSMXaSIOESLZP ANGLE-MAX 
7-SSMHaSIOESLIP AHGLE-MXH 
1-SSSOaS I OESLI P ANGLE— SO 
9-NTA •NACELLE TILT ANGLE 

0- PHOTaCLIMB/OESCENT HATE 

1- PCAS-CABIN ?HES SURE-MEAN 

2- PCMX-CABIN PHESSUHE-MAX 

3- PCMHaCABIH PHESSUHE-HXN 

4- PCSQaCABIH PHESSUHE-SO 
9-COP -CABIN 01 PP PHESSUHE 

5- TCABaCABIH TEMPEHATVHE 
7-CC •CABIN SPEED OP SOUHO 
•-SIGC-CABIN DENSITY HATZO 
•—PC -CABIN AIH DENSITY 
O-PMCC-CABIN CHAH . IMPEDANCE 


4 

XO 

xo 

XO 

xo 


TS 

xo 

10 

xo 

xo 

xo 

XX 

XX 

XX 

XX 

XI 


X-TQEi-POWEH TUHBINE TOHqUE-MEAH 

2- TQMX-ORXVE ENG PWR TUH8 TORQUE-MAX 

3- TQKN-ORXVE ENG PWR TUHB TORQUE-MXN 


-TOSOaOHIVE ENG PWR TUHB TORQUE-ID 
- PSHPa POWEH TUHBINE SHP 
-SHPC-OHIVE ENG POWEH TUHBINE %MCP 

-HPNL-POWEH TUHBINE % HP-MEAN 

-HLMXaOHZVE ENG PWH TUHBINE % HP- MAX 

-HLMHaORIVE ENG PWH TUHBINE % HP-MXH 

-HLSOaOHZVE ENG PWH TUHBINE % NP-SD 

-PRCfapWH SECTION MAX CONTINUOUS PWH 
-HPT aOHZVE ENG PWH TUHB ROT SPO-HP 
-RPNNaCOMPHESSOH % NO-MEAN 
—RHMX- DRIVE EHG COMPRESSOR % NG-MAX 
-HHMHaOHZVE EHG COMPRESSOR %.NG-MXN 
-RNSOaOHIVE ENG COMPRESSOR % NG-SO 
RPNC-OXIVE EHG COMP ROTATION SPO-NG 
-CPX aOHIVE EHG COMP XST ROTOR PREQ 
-Tft -OHIVE ENG TUHB LAST ROTOR fRCQ 


RPM3 

-R3MX* 

-R3MN 

-R3S0* 

-RPMX- 

— RXMX 

— RXMH 
-RXSD 
-SPHXa 
-L SPX« 
— LST1- 


LP COMPRESSOR 
LPT SPEY EHG 
LPT SPEY ENG 
LPT SPEY EHG 
HP COMPRESSOR 
LPT SPEY EHG 
LPT SPEY ENG 
LPT SPEY ENG 
LPT SPEY ENG 
LPT SPEY ENG 
LPT SPEY ENG 


Ill-RPM4«LP COMPRESSOR SPO-NX 4 
Xli-R4MX»RT SPEY ENG LP COMP SPEED-MAX 
XX7-R4MNaRT SPEY ENG LP COMP SFCEO-MXN 
X1B-JI4JO-RT SPEY ENG LP COMP SPEEO-SO 
XXS-RPM2-MP COMPRESSOR SPO-N2 % 
XZ0-R2MX-RT SPEY ENG HP COMP SPEED-MAX 
XZX-RlMHaRT SPEY ENG HP COMP SPCEO-NIN 
122- RISO-RT 5 PE7 ENG HP COMP SPEED-SO 
X23-SPn:-PIGHT SPEY CMC NET THRUST 
124-R3FX-P.T SPEY ENG LPC:XST ROT BPF 
X 2 3-RST X •RT SPEY ENG LPT: LAST ROT BP P 


<D 

C 

•H 

c 

Ci3 

Cse 

CL CD 
> 
•H 
U 

a 


-PRPMa PROP PAH ROTATION SPEED 
—VROT* PROPFAH TIP ROTATION VEL 
-VFWOaPROPFH TIP fWO VEL-TRU AIH SPO 
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A min or modification was made to the SRPIE (Single Rotation Propfan Installation 
Effects) code to enable the calculation to be performed at specified values of 4> (the azimuthal 
angle), rather than at predetermined 15° intervals. 

No changes were made to the SRPFS (Single Rotation Propeller Fuselage Scattering) 
code. 

3.3 Results 

3.3.1 Case 1: Design Point 

3.3.1. 1 Uninstalled Predictions 

Results calculated under free field conditions at the fuselage microphone locations along 
the 93.65 and 91.71 water lines and at the boom microphones are shown in Figures 81 and 82, 
respectively. The first three harmonics of BPF are shown in each figure. 

It will be noted that, under these conditions, the predicted levels are the same at the 
fuselage surface (Figure 81) and at the boom (Figure 82). This is what would be expected from 
the prediction model, given that the boom and fuselage microphones were at virtually the same 
sideline distance (Table 5), and that the model predicts no azimuthal variations. 

One feature depicted in Figure 81 (which is missing from Figure 82, due to the absence of 
amicrophone at that location), is that all three harmonics of BPF are predicted to peak forward 
of the plane of rotation. In order to examine this feature, Figure 83 was plotted, showing the 
BPF tone predicted by the steady loading and thickness components of the acoustic model 
separately, together with the total already seen in Figure 81. It appears (from Figure 83) that 
the thickness noise prediction, combined with phasing effects, is causing the peak in the 
forward arc, with the relative phasing between the two components actually reducing the total 
level below the peak predicted by the steady loading model at Fuselage Station 355. 

One final effect, seen in Figure 84, concerns the selection of the angle used to divide the 
lift on any blade section into its thrust and torque components. As described in Section 2.1, 
there is an option in the SRPFAN code allowing’ the use of either the helicoidal angle (p A , in 
Figure 1) or the blade pitch angle (Pn> i n the same figure) in the evaluation of ky in Equa- 
tion 2.1.2. Figure 84 indicates the minor differences found between the two models. In 
general, use of the helicoidal surface option results in levels that are 1 to 2 dB higher than 
those obtained with the blade pitch angles. All of the results presented in this report (with the 
exception of Figures 84, 96, and 100), use the blade pitch angles to divide the section lift into 
its thrust and torque components. 

3.3.1 .2 Installation Effects 

The effects of the installation environment (presence of a fuselage and wing lifting line, 
together with angle-of-attack operation) were calculated, again under free field conditions, at 
the microphones located on the left side of the fuselage in the plane of the propfan (Fuselage 
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Figure 81. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 1) 
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Figure 82. PTA Airplane Boom Predictions - Uninstalled (Prediction Case No. 1). 
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Figure 83. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 1). 
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Figure 84. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. I). 





Station 301) and 0.25 propfan diameters upstream and downstream of this plane, at Fuselage 
Stations 274 and 328, respectively. Results are listed in Table 10 for the first and second 
harmonics of BPF, and the free field installed levels of the BPF tone are shown graphically in 
Figure 85. The effect of the installation is predicted to be small everywhere. 

In order to make the prediction, it was assumed that the inflow angle to the propfan was 
given by the parameter “LA” (the combination of the airplane angle-of-attack and the nacelle 
tilt) in Table 6. For this case the value is low, namely 0.81°. Other inputs required are the 
wing lift coefficient and the location of the center of the propfan, relative to the wing lifting 
line. The lift coefficient was calculated from the aircraft gross weight and other relevant 
parameters given in Table 6, together with an estimate of the wing area scaled from Figure 78. 
In addition, the analytical model employed in the SRPIE code for the flowfield of a wing lift- 
ing line is restricted to rectangular (unswept) wings. 

Consequently, it was assumed that the unswept wing lifting line passes through the 
midchord of the wing/fuselage installation. Finally, the blade section lift and drag coefficients 
versus angle-of-attack data required by the prediction code were obtained using the perform- 
ance data given in Reference 33. 1 for the SR-7A model propfan in the CLDSRP performance 
preprocessor. 

3.3.1 .3 Fuselage Scattering Effects 

The final set of calculations performed for this high flight Mach number case were con- 
cerned with the effects of refraction and diffraction described in Section 23. These effects 
were shown to increase with increasing flight Mach No.; hence it was decided that, for this 
case only, an assessment would be made of their contribution. 

Results were calculated for three values of fuselage boundary layer thickness: 8 = 1 inch, 
8 = 3 inches (suggested by Lockheed), and 8 = 6 inches. The predicted effects are shown in 
terms of “dBFW,” which is defined as: 


dBFW = 201og 10 (P(Rf,4>)/Pi(Rf+8,4>) 

where P(R,4>) = Pi(R,<|>) + Ps(R,4>), (incident + scattered fields.) 

For pure reflection, dBFW = 6 dB. 

Figure 86 presents the predicted effects of the three different boundary layer thicknesses 
on the fundamental BPF tone as a function of distance along the fuselage at a fuselage 
azimuthal angle, 4 >-fuse (Table 5), equal to zero. This approximately corresponds to the 93.65 
and 91.71 water line locations used in the uninstalled predictions. The results are similar to 
those presented in Section 23, with the greatest decreases relative to the pure reflection 6 dB 
value occurring ahead of the plane of the propfan (Fuselage Station 301) as a function of 
boundary layer thickness, and the boundary layer being predicted to have negligible effect 
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Table 10 . Predicted Installation Effects for Case 1 (Free Field). 
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Figure 85. PTA Airplane Fuselage Prediction - With Installation Effect (Prediction Case No. 1) 
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Figure 86. PTA Airplane Fuselage Scattering Prediction (Case No. 1). 


downstream of approximately 0.25 propfan diameters aft of the plane of rotation (Fuselage 
Station 328). 

Figure 87 shows the result of combining these effects with the BPF tone of Figure 81. It 
now appears that the location of the peak level is a function of fuselage boundary layer thick- 
ness, with the high level in the forward arc diminishin g as the thickness is increased, while the 
levels aft of the propfan remain constant 

Figures 88 through 90 illustrate the effects predicted on the BPF tone at points around the 
left side of the fuselage, of the same three values of fuselage boundary layer thickness (1 inch, 
3 inch, and 6 inch) as before, at Fuselage Stations 274, 301, and 328, respectively, in terms of 
dBFW. The “y-coordinate” is “^fuselage,” where zero represents the position on the fuselage 
of a line connecting the midpoint of the propfan disk to the centerline of the fuselage. 

Finally, Figures 91 through 93 represent the results of combining the fuselage scattering 
effects shown in Figures 88 through 90 with the installation effects of Table 10 and the isolated 
propfan steady loading and thickness prediction, to give predictions of the BPF tone levels 
that would be measured at the microphone locations around the fuselage in the plane of rota- 
tion of the propfan (Figure 92), and 0.25-diameters forward and aft of that plane (Figures 91 
and 93, respectively). 

These figures demonstrate the effects of the increase in fuselage boundary layer thickness, 
both along and around the fuselage surface. For example, in Figure 93 it can be seen that, 
whereas at the measurement locations at and below the 91.71 water line the predicted tone 
level is virtually independent of boundary layer thickness, above this point there is a relation- 
ship. It will be remembered (Figure 77) that the propfan rotates clockwise when viewed from 
forward looking aft The blades are, thus, approaching the fuselage below the 91.71 water line, 
and retreating above it Both the installation effects (other than angle-of-attack) and the 
fuselage scattering results are functions of the direction of rotation. 

3.3.2 Case 2: Low Altitude, Full Power, and Tip Speed 
3.3.2. 1 Uninstalled Predictions 

The input parameters used for these predictions are listed in Table 7, and the results 
predicted for the first three harmonics of BPF under free field conditions along the 91.71 and 
93.65 water lines on the fuselage, and at the boom microphone locations, are shown in Figures 
94 and 95. Comparison of these figures with Figures 81 and 82, which show equivalent predic- 
tions for the design point of the propfan, demonstrates the different character of the predicted 

tone noise from a propfan at supersonic helical Mach number (Case 1: Mht * 1.137), and 

from the same propfan at a helical tip Mach No. of 0.787. At the lower Mach number, the 
directivity patterns for the three harmonics are smoother, with the predicted peak level occur- 
ring behind, or at, the plane of rotation. The peak predicted level is 10 to 12 dB lower at the 
lower tip helical Mach number, and the decrease in tone level with increase in harmonic order 
is more pronounced. The harmonic levels of Figure 81 are probably held up by the thickness 
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Figure 87. PTA Airplane Fuselage Prediction - With Fuselage Scattering (Prediction Case No. 1). 




100 


O O 


asnj - i.iid 


142 


O 

m 

i 


00T- 






Delta 


n 

a 

u 





a 

■*4 

m 

Q 



esn j - md 


o 

CN 



c 

o 

2 

oo 


1 


<D 

O) 

0) 

CO 

D 

Li- 


03 


O J^r 
CM CD 
I "O 


O 

00 


o 

I 


143 


Figure 89. PTA Airplane Fuselage Scattering Prediction (Case No. I) 
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Figure 91. PTA Airplane Fuselage Prediction - with Fuselage Scattering (Prediction Case No. I 
Fuselage Station 274). 
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Figure 92. PTA Airplane Fuselage 
Fuselage Station 301). 
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Figure 93. PTA Airplane Fuselage Prediction - With Fuselage Scattering (Prediction Case No. 1 
Fuselage Station 328). 
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Figure 95. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 2). 





noise prediction which, as seen in Figure 83, is of equal or greater importance at the high helical 
tip Mach number point In Figure 94, the decrease in level with increase in harmonic number 
is consistent with noise due to steady loading as the dominant source. As before, free field 
predictions at the boom microphone locations (Figure 95) are virtually identical to those illus- 
trated in Figure 94. 

As in the previous case, the two options for splitting the lift on the blade sections into their 
thrust and torque components were investigated. The results are shown in Figure 96. It is 
interesting to note, in comparing Figure 96 with Figure 84, that the two methods have “switched 
positions." For this low Mach number, off-design case, the predictions using the helicoidal 
surface in the acoustic prediction are consistently 1 to 2 dB lower in the peak region than those 
obtained with the blade pitch angles. 

3.3^2 Installation Effects 

Installatio n effects for this case were calculated using an angle -of-attack of 3.7°, obtained 
from Table 7. The wing lift coefficient was calculated from the aircraft gross weight and other 
relevant parameters from the same table. No low flight Mach number performance data was 
available for the SR-7 propfan design, so it was assumed that the Q, and C D versus angle-of- 
attack characteristics of the blade sections were s imil ar to those of the SR-3 model propfan. 

Accordingly, data from the performance maps of Reference 33.2 were used as input to 
program CLDSRP. The results of the calculation are tabulated in Table 11 for the first and 
second harmonics of BPF, and Figure 97 shows the free field BPF tone levels predicted around 
the left side of the fuselage with installation effects included in the prediction. 

Once again, the predicted installation effects are small, with the predicted decrease toward 
the top of the fuselage and increase toward the bottom of the fuselage in line with what would 
be anticipated with a positive angle-of-attack. No fuselage scattering predictions were made 
for this case. 

3.3.3 Case 3: Low Altitude, Reduced Power, and Tip Speed 
3.3.3.1 Uninstalled Predictions 

Figures 98 through 100 are equivalent to Figures 94 through 96 for the full power and tip 
speed case described above. Examination of Tables 7 and 8 shows the predicted difference in 
the free field peak levels of Figures 94 and 98 to be a result that relates more to the difference 
in tip Mach number than to the difference in power absorbed in the two configurations. 

The tip speed effect on steady loading noise from a point force can be written: 

A(dB) = 20*Log 10 (V2/Vl) ~ 6 
= 7.9 dB in this case; 
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Figure 96. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 2) 





Table 11. Predicted Installation Effects for Case 2 (Free Field). 
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Figure 97. FTA Airplane Fuselage Prediction - Willi Installation EiTects (Prediction Case No. 2) 
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Figure 98. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 3). 
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Figure 100. PTA Airplane Fuselage Predictions - Uninstalled (Prediction Case No. 3). 


while, for the difference in power: 

A(dB) * 20*Logio(shp2/shpl) 

= 1.9 dB. 

The total of 9.8 dB is very close to the 9.4 dB predicted by the distributed source model, 
implying that, under these low speed conditions, the added complexity of the full model is not 
necessarily needed. 

Once again, the boom results of Figure 99 are very similar to those calculated at the 
microphone locations on the fuselage (Figure 98). For this case, the nacelle tilt was +2°, 
leading to a slightly greater difference in sideline distance between the boom and the fuselage 
than was seen in the previous two cases (Table S). 

The differences between results calculated using the helicoidal surface angles to divide 
the lift into its components for the acoustic calculation and those obtained with the blade pitch 
angles are portrayed in Figure 100. As in Figure 96, the helicoidal surface results are approxi- 
mately 2 dB lower than those obtained with the blade pitch angles. 

3.3.3.2 Installation Effects 

The nacelle tilt of +2° for this case resulted in an angle-of-attack of 5.94° (Table 8). The 
other installation inputs were calculated as before, from parameters also in Table 8. 

The free field predicted results are shown in Table 12 and Figure 101. Although slightly 
larger than those predicted for Case 2, the predicted effects are again small, with the increase 
toward the bottom of the fuselage and the decrease toward the top in accordance with the 
expected angle-of-attack effect Overall, the predicted free field levels around the fuselage 
(Figure 101) are lower than those shown for Case 2 (Figure 97) and exhibit less variation with 
fuselage station. Again, no fuselage scattering effects were calculated for this case. 

3.4 Conclusions 

The acoustic models developed under this contract as described in Section 2.0 of this report 
for the prediction of high speed turboprop noise resulting from steady loading and thickness 
sources on the blading together with installation effects arising from: angle-of-attack opera- 
tion, the presence of a cylindrical fuselage, and the presence of a wing lifting line, have been 
exercised to predict the free field noise at various microphone locations on the aircraft for 
three flights of the Propfan Test Assessment airplane. In addition, for one flight representing 
the design point of the propfan, the fuselage scattering code which was also developed under 
this contract and described in Section 2.0 of this report was utilized to predict the levels that 
would be measured by microphones on the fuselage of the airplane, given any of three pos- 
sible fuselage boundary layer thicknesses. 
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Tabic 12. Predicted Installation Effects for Case 3 (Free Field). 
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Figure 101. PTA Airplane Fuselage Prediction - With Installation Effect (Prediction Case No. 3) 


The results show that at high' helical tip Mach number (Case 1), the noise predicted due 
to the thickness (or volume) of the blades is comparable to, or even greater than, that predicted 
to result from the steady loading component alone. The slow dropping off of the peak levels 
with increase in harmonic number for this case has been ascribed to the thickness source. 

The inclusion of fuselage scattering effects demonstrates that the location of the peak level 
on the fuselage at cruise can be a function of the thickness of the boundary layer on the fuselage. 
Under free field conditions, the maximum level is predicted to occur forward of the propfan 
plane. It has been shown that the presence of the boundary layer leads to a predicted reduction 
in the tone levels on the fuselage upstream of the propfan, while giving almost pure reflection 
in the aft region. 

Installation effects at anise are predicted to be less than 1 dB over the region in the vicinity 
of the propfan plane. 

In the low flight Mach number regime (Cases 2 and 3) the predicted noise levels are 
dominated by steady loading effects. The predicted differences between Cases 2 and 3 (“full 
power and tip speed” versus “reduced power and tip speed”) are shown to result more from 
the tip speed effect than from the reduction in power absorbed. 

The predicted effects of the installation indicates levels increasing toward the bottom of 
the fuselage and decreasing toward the top in both cases (which is consistent with the angles- 
of-attack input), although the levels of the predicted effect are low everywhere. 

An examination of the consequences of the use of the blade pitch angles rather than the 
helicoidal surface angles in dividing the blade section lift into its thrust and torque components 
for the acoustic calculation shows the difference between the two sets of results to be of the 
order of 1 to 2 dB. One interesting point noticed is that: whereas at cruise, results generated 
using the helicoidal surface are higher, in the low flight Mach number regime, the opposite is 
true. It is not known if this is a universal result, or a function of the particular geometry under 
consideration here. 


160 



4.0 OVERALL CONCLUSIONS 


Acoustic models have been formulated and programmed under this contract to predict 
the tone noise generated by a high speed turboprop as a result of: the steady loading on the 
blades; the thickness (volume) of the blades; and the unsteady blade loads resulting from the 
installation environment A model to predict modifications to the sound field emitted from 
the turboprop resulting from the curvature of an airplane fuselage and its attached boundary 
layer also has been formulated and programmed. 

The steady loading and thickness model matches a frequency-domain, noncompact source, 
linear acoustic model with the solutions from a nonlinear, aerodynamic flow field analysis 
program that are used to obtain the blade loading distributions required as input by the acoustic 
code. The model was developed originally for the acoustic far field, and a semiempirical 
adjustment has been included to account for near field effects. An alternate method for 
predicting near field noise directly from the numerical flow field solution was tried and 
discarded as impractical at this time. 

Results from the steady loading and thickness model have been compared with data from 
three model propfans (SR-2, SR-3, and SR-6) operating at high flight speed, both in a wind 
tunnel environment and in flight (mounted on the Jetstar aircraft). In general, the predictions 
show good agreement with the Jetstar boom data and moderate agreement with the data 
measured on the airplane fuselage and in the wind tunnel. 

The installation effects model developed under this contract employs a quasi-steady 
formulation to predict the thrust and tangential forces on the blades resulting from the non- 
uniform flow environment in which the propfan operates when installed on an aircraft. The 
specific installation effects considered here are: angle-of-attack operation, the presence of a 
fuselage, and the presence of a wing lifting line. Calculation of the flow fields resulting from 
these effects has been included in this work. 

The acoustic model is less sophisticated than that employed for the steady loading and 
thickness calculation, being compact in the chordwise direction, so the results obtained are 
presented as a delta, which can be applied to the outcome of a steady loading and thickness 
calculation. 

Comparisons with SR-2 model data taken at a low wind tunnel Mach number show encour- 
aging agreement. 

The final model developed under this contract concerns the scattering of the sound field 
emitted by the propfan as a result of the presence of an airplane fuselage with its associated 
boundary layer. In the analysis, the emitted sound field is expressed in terms of multiple 
cylindrical waves converging onto the cylindrical fuselage. These waves are then scattered 
from the surface of the cylinder, or if a boundary layer is present, the effects of refraction 
through this layer are included. The results obtained with this model are in agreement with 
observed trends. 
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The three models developed under this contract have been exercised for three flights of 
the Propfan Test Assessment airplane. In this project, a 9-foot-diameter SR-7L propfan was 
mounted on the left wing of a modified Gulfs tream n business jet Microphones were installed 
along and around the fuselage, as well as on a wing-mounted boom. The predictions were 
made for flights whose operating conditions corresponded to the design point (high altitude 
cruise); a low altitude, high power and tip speed point; and a low altitude, reduced power and 
tip speed point Data comparisons were not performed, but free field levels were predicted 
at various microphone locations for all three cases. 

Results indicated the thickness component to be a strong contributor at cruise, with the 
location of the maximum level on the fuselage a function of the thickness of the fuselage bound- 
ary layer. In the low altitude, low flight Mach number regime, steady loading sources were the 
major contributors, with tip speed effects dominant Installation effects were predicted to be 
fairly small in all three cases, with angle-of-attack effects appearing to dominate. 

In conclusion therefore, three major computer programs have been delivered under this 
contract They predict the noise generated by a high speed propfan as a result of its steady 
loading and thickness and its installation environment together with corrections that should 
be applied to measurements taken on the surface of a fuselage. 

These programs have been compared with data from a variety of sources, and on the whole, 
demonstrate reasonable agreement In addition, two smaller programs have been delivered; 
one to convert the results of an aerodynamic flow solution to the input format required by the 
steady loading and thickness prediction code; the other to generate blade section Cl and Cd 
versus angle-of-attack data for input to the installation effects code. Possibilities for further 
work include the development of a rigorous near field model, and the effect (in the far field), 
of the presence of a fuselage in the near field. 


i 
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